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STUDIES  IN  TURBULENCE  AND  TURBULENCE  CONTROL 


The  following  projects  have  been  completed.  A  brief  discussion  of  the  work  is  given,  and  the 
principal  publications  are  listed. 

A.  Projects  which  bear  on  research  in  reacting  flows 

1.  Differential  diffusion:  Consider  two  scalars  homogeneously  mixed  in  the  plenum  chamber 
of  a  jet.  The  jet  is  turbulent.  The  two  scalars  have  different  Schmidt  numbers  (i.e.,  their  molecular 
diffusion  coefficient  are  different).  What  effect  will  this  difference  in  Schmidt  numbers  have  in  the 
far  field  in  the  instantaneous  spatial  and  temporal  distributions  of  the  two  scalar  fields?  This  is  a 
fundamental  issue  in  the  context  of  combustion  where  several  species  with  different  diffusion 
properties  are  mixed  by  turbulence.  A  popular  assumption  made  in  combustion  studies  is  that 
differential  diffusion  effects  are  negligible. 


Experimental  data  on  differential  diffusion  between  two  species  with  large  and  quite  disparate 
Schmidt  numbers  were  obtained  in  a  turbulent  water  jet  by  optically  measuring  the  two  species 
concentrations  simultaneously.  Experimental  conditions  were  chosen  to  eliminate  possible  inertial 
effects.  Schmidt  number  ratios  of  4  and  18  were  considered.  Differential  diffusion  was  found  to  be 
statistically  significant  and  to  manifest  at  scales  far  larger  than  the  Batchelor  scale.  In  some 
instances,  the  concentration  signal  for  the  species  with  larger  diffusivity  was  simply  a  blurred 
version  of  the  other,  while  in  some  instances  structures  present  in  one  signal  were  absent  from  the 
other.  This  second  observation,  presumably  a  more  complex  effect  due  to  diffusion  across  velocity 
gradients,  has  been  discussed. 

The  work  is  submitted  for  publication  to  the  Journal  of  Fluid  Mechanics  as:  Differential 
diffusion  in  low  Reynolds  number  water  jets,  by  J.R.  Saylor  and  K.R.  Sreenivasan 


.13 


I, 

ir'  '  *  - 


.  or 


rrj 


3.  P  hotobleaching  of  sodium  fluorescein  in  water:  As  a  secondary  outcome  of  the  differential 
diffusion  study,  data  concerning  the  photobleaching  of  sodium  fluorescence  in  water  were 
obtained.  The  importance  of  avoiding  this  effect  in  fluorescence  studies  is  emphasized. 

The  results  are  available  in  preprint  form:  The  temporal  behavior  of  the  photobleaching  of 
disodium  fluorescein  in  water,  by  J.R.  Saylor 

2.  Scalar  dissipation  conditioned  on  the  scalar  concentration:  In  the  pdf  methods  for  the 
scalars,  the  quantity  that  needs  modeling  is  the  conditional  expectation  of  the  scalar  dissipation 
conditioned  on  the  scalar  concentration  itself.  To  date,  there  are  no  reliable  measurements  of  this 
quantity.  We  have  measured  it  in  the  following  cases:  dye  concentration  fluctuations  in  a  water  jet, 
temperature  fluctuations  in  the  wake  behind  a  heated  cylinder,  and  temperature  fluctuations  in  high- 
Reynolds-number  atmospheric  turbulence.  The  effects  of  Taylor's  hypothesis  are  tested  by 
obtaining  both  spatial  and  temporal  data.  Similarly,  the  effect  of  probe  resolution  has  also  been 
tested  by  making  measurements  at  different  resolutions.  Finally,  the  effect  of  using  one  component 
of  the  scalar  dissipation  instead  of  the  sum  of  all  its  three  components  has  been  studied. 

This  work  has  been  has  been  submitted  for  publication  to  Physics  of  Fluids  A:  Conditional 
scalar  dissipation  rates  in  turbulent  wakes,  jets  and  boundary  layers,  by  P.  Kailasnath,  K.R. 
Sreenivasan  and  J.  Saylor. 

3.  Helium  jets:  We  have  finally  completed  a  detailed  study  of  with  helium  jets  issuing  into  air. 
'The  principal  results  were  summarized  in  the  final  report  on  AFOSR-87-Ol  16,  and  the  publication 
summarizing  the  results  has  appeared  in  the  Journal  of  Fluid  Mechanics  (vol.  249,  pp.  619-664, 
1993),  and  a  copy  is  enclosed. 

4.  The  thickness  distribution  of  the  OH  regions  in  turbulent  diffusion  flames:  The  OH 
radicals  in  a  hydrogen  diffusion  flame  burning  in  ambient  air  were  visualized  using  planar  laser- 


induced  fluorescence.  Local  thickness  of  the  OH  regions  were  measured  at  three  downstream 
distances.  Measurements  show  that  the  thickness  is  distributed  approximately  lognomnaUy. 

The  results  were  published  as:  The  thickness  distribution  of  Of!  regions  in  a  turbulent  diffusion 
flame,  by  A.  Johnson,  K.R.  Sreenivasan  and  M.  Winter,  Comust.  Sci.  and  Tech.  89,  1  4, 
1992.  A  copy  is  enclosed. 

B.  Basic  turbulence  work 

1.  Kolmogorov's  refined  hypotheses:  In  1962,  Kolmogorov  put  forward  a  vital  refinement 

of  his  earlier  phenomenological  theory  of  high-Reynolds-number  turbulence.  An  important 
quantity  in  this  description  is  the  flux  of  energy  <j>r  transferred  across  scales  of  size  r.  Kolmogorov 

assumed  that,  in  the  inertial  range,  <j)r  is  the  only  relevant  quantity  upon  which  all  other  quantities 
would  depend.  Furthermore,  he  identified  4>r  with  rer,  where  Er  is  the  rate  of  energy  dissipation 

per  unit  mass  averaged  over  a  volume  of  linear  scale  r.  The  basis  of  Kolmogorov's  theory  is 
contained  in  the  following  two  celebrated  hypotheses. 

The  first  similarity  hypothesis:  If  r«L,  where  L  is  a  measure  of  the  large  scale  of  turbulence, 
the  probability  density  function  (pdf)  of  the  stochastic  variable 

v  Au(r) 

(rer)1/3 

depends  only  on  the  local  Reynolds  number  Rer  -  r(rer)1/3/V,  where  v  is  the  kinematic  viscosity 

of  the  fluid  and  Au(r)=u(x+r)-u(x),  u  being  the  x-component  of  the  velocity  vector  u(x)  and  r  is 
measured  along  x. 

The  second  similarity  hypothesis:  If  Rer»l,  the  pdf  of  V  does  not  depend  on  Rer  either  (nor 
on  r,  and  is  therefore  universal). 

We  have  earlier  shown  that  the  pdf  of  V  depends  on  r  as  well,  but  also  showed  that  several 

other  aspects  of  these  hypothecs  hold.  F.vcn  before,  this  verification,  consequences  of  these 


hypotheses  had  been  used  extensively  in  the  turbulence  literature.  In  spite  of  their  widespread  use, 
the  hypotheses  pose  some  troubling  problems;  the  primary  shortcoming  is  that  they  have  not  yet 
been  derived  from  basic  principles.  Under  the  assumption  of  statistical  independence  of  velocity 
increments  across  scales  of  the  order  of  the  Kolmogorov  scale,  we  have  been  able  to  show  that  a 
modified  version  of  Kolmogorov's  refined  similarity  hypotheses  follows  purely  from  probabilistic 
arguments.  The  connection  of  this  result  to  three-dimensional  fluid  turbulence  has  been  discussed. 
It  has  been  shown  that  the  physical  picture  of  a  cascade  need  not  be  assumed  a  priori ,  but  rather  as 
an  a  posteriori  interpretation. 

The  work  will  appear  in  an  AIP  book  as:  Independent  velocity  increments  and  Kolmogorov's 
refined  similarity  hypothesis,  by  G.  Stolovitzky  and  K.R.  Sreenivasan 

2.  Self-affine  time  traces  of  turbulent  signals:  Consider  a  snapshot  of  a  turbulent  flow,  and 
define  iso-velocity  or  iso-concentration  surfaces  in  that  snapshot.  Previous  studies  have  shown  that 
such  surfaces  are  self-similar  and  can  be  characterized  by  fractal  dimensions.  In  contrast,  consider 
a  time  trace  of  temperature  or  velocity  at  a  single  point  in  a  turbulent  flow.  For  such  quantities,  the 
abscissae  (i.e.,  time)  and  the  ordinate  (i.e.,  temperature  or  velocity)  are  two  very  dissimilar 
entities,  and  the  signal  can  be  rendered  into  a  variety  of  objects  by  stretching  one  or  the  other  axis. 
For  example,  by  stretching  the  time  axis,  one  can  make  the  signal  look  very  smooth;  it  will  then 
have  a  dimension  of  unity.  By  compressing  the  time  axis,  the  signal  can  be  made  very  rough  and 
almost  space-filling.  These  curves  are  therefore  harder  to  characterize  in  terms  of  fractal 
dimensions.  We  have  finally  understood  the  way  in  which  this  can  be  done. 

The  material  is  written  up  for  publication  as:  Fractal  dimensions  of  time  series  in  turbulent 
flows,  by  K.R.  Sreenivasan  and  A.  Juneja,  submitted  to  J.  Fluid  Mech.,  1992 

3.  Synthetic  turbulence:  On  the  basis  of  the  fractal  and  multifractal  work  with  which  we  have 
been  concerned  so  far,  we  have  been  able  to  outline  a  scheme  for  constructing  a  stochastic  field  that 
differs  in  no  significant  way  statistically  from  velocity  fluctuations  in  high-Reynolds-number 


turbulence.  We  designate  this  field  as  synthetic  turbulence.  The  scheme  for  generating  synthetic 
turbulence  consists  of  the  following  three  ingredients,  implemented  here  in  one  dimension.  From 
such  signals,  one  can  construct  velocity  increments,  box-averaged  dissipation  rate,  Kolmogorov's 
universal  variable  V,  all  of  which  have  the  same  properties  as  their  counterparts  in  the  real  signal. 
Also,  properties  such  as  the  power  spectral  density,  correlation  function,  fractal  dimensions,  etc 
can  be  duplicated.  The  importance  of  this  finding  is  two-fold: 

(a)  It  shows  that  the  physical  picture  we  have  of  the  turbulence  dynamics  in  high-Reynolds- 
number  turbulence  is  essentially  correct  up  to  some  level  of  accuracy. 

(b)  If  one  can  generate  a  field  that  is  so  close  to  turbulence,  such  a  field  can  serve  as  good 
initial  conditions  for  direct  numerical  simulations. 

The  work  is  available  in  the  form  of  a  preprint  as:  Synthetic  turbulence,  by  A.  Juneja,  D.P. 
Lathrop,  K.R.  Sreenivasan  and  G.  Stolovitzky. 


4.  The  wavelet  decomposition  of  turbulence:  The  wide  range  of  spatial  and  temporal  scales 
excited  in  fully  turbulent  flows  have  been  historically  represented  by  Fourier  modes.  It  has  been 
recognized  recently  that  wavelets  possess  several  advantages  over  Fourier  modes,  especially  in 
strongly  inhomogeneous  shear  flows  dominated  by  localized  structures.  We  apply  wavelet 
transform  to  independent  realizations  of  two-dimensional  sections  of  the  dye  concentration  field 
in  turbulent  water  jets  at  moderate  Reynolds  numbers,  and  show  that  the  smaller  the  scale  the 
more  stringy  it  is  in  two-dimensional  intersections.  The  transformation  from  the  large  scale  to 
small  scale  can  be  better  described  by  a  self-affine  process  associated  with  the  stretching  due  to 
the  velocity  field  rather  than  by  a  striedy  self-similar  process  characterized  by  a  geometrically 
similar  break  down  process.  Statistical  measures  such  as  the  number  of  structures  at  a  given  scale 
and  the  average  aspect  ratio  of  these  structures  are  calculated.  We  also  analyze  the  space-time 
data  on  two-dimensional  sections  of  the  concentration  field,  and  analyze  interactions  of  various 
wavelet  scales.  A  principal  conclusion  to  emerge  from  the  analysis  is  that  the  wavelet 
coefficients  which  are  disparate  in  scale  or  separated  widely  in  space  are  nearly  statistically 
independent.  This  notion  has  been  put  to  use  in  predicting  the  probability  density  function  of 
wavelet  coefficients  (wpdf).  A  theory  has  been  worked  out  for  the  wpdf  for  passive  scalars  in 
tuihulcnt  flows,  and  the  resulting  equations  arc  closed  by  invoking  the  Fokkcr-Planck 


approximation  as  well  as  the  statistical  independence  of  wavelet  coefficients.  The  theory  is 
compared  with  the  results  from  experiments.  It  is  found  that  the  theory  works  well  for  large  and 
intermediate  scales  but  not  for  the  small  scales.  The  reasons  for  this  failure  of  the  theory  have 
been  explored,  and  it  is  shown  that  a  particular  form  of  statistical  dependence  makes  the  theory 
in  agreement  with  the  experiment. 

One  paper  has  been  written  and  submitted  to  Physics  of  Fluids  A:  The  wavelet  coefficient 
probability  density  function  for  turbulent  flows,  by  P.L.  Similon  and  K.R.  Sreenivasan.  A 
seconfd  paper  is  available  in  draft  form  as:  The  multiscale  structure  of  turbulent  water  jets,  by 
K.R.  Sreenivasan  and  R.M.  Everson 

5.  An  update  on  the  intermittency  exponent  in  turbulence:  The  issue  of  the  experimental 
determination  of  the  intermittency  exponent  ,  m,  is  revisited  and  it  is  shown  that  the  ’best’ 
estimate  for  it  is  0.2510.05.  This  ’best’  estimate  is  obtained  from  recent  atmospheric  data,  and  is 
based  on  several  different  techniques  of  measurement.  The  work  appeared  in  Phys.  Fluids  A,  and 
a  copy  is  enclosed. 

6.  Other  work:  Aspects  such  as  the  scaling  of  the  structure  functions  in  the  inertial  as  well  as 
dissipative  scale  ranges,  the  effect  of  filtering  on  the  Kolmogorov’s  4/5  relation,  the  scaling 
relations  of  circulation,  have  been  studied  in  some  detail.  They  will  not  be  reported  here.  As  and 
when  the  results  appear  in  print,  copies  of  the  papers  will  be  mailed  to  AFOSR. 
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A  study  has  been  made  of  the  instability  and  the  subsequent  breakdown  of 
axisymmetric  jets  of  helium/air  mixtures  emerging  into  ambient  air.  Although  the 
density  of  the  nozzle  gas  is  less  than  that  of  the  ambient  fluid,  the  jet  is  essentially  non- 
buoyant.  Two  kinds  of  instability  are  observed  in  the  near  field,  depending  upon  the 
mean  flow  parameters.  When  the  ratio  of  the  exiting  nozzle  fluid  density  to  ambient 
fluid  density  is  f\J/> ,  >  0.6.  shear-layer  fluctuations  evolve  in  a  fashion  similar  to  that 
observed  in  constant-density  jets:  the  power  spectrum  near  the  nozzle  is  determined  by 
weak  background  disturbances  whose  subsequent  spatial  amplification  agrees  closely 
with  the  spatial  stability  theory.  When  the  density  ratio  is  less  than  0.6.  an  intense 
oscillatory  instability  may  also  arise.  The  overall  behaviour  of  this  latter  mode  (to  be 
called  the  'oscillating'  mode)  is  shown  to  depend  solely  upon  the  density  ratio  and 
upon  D/H,  where  D  is  the  nozzle  diameter  and  0  is  the  momentum  thickness  of  the 
boundary  layer  at  the  nozzle  exit.  The  behaviour  of  this  mode  is  found  to  be 
independent  of  the  Reynolds  number,  within  the  range  covered  by  the  present 
experiments.  This  is  even  true  in  the  immediate  vicinity  of  the  nozzle  w  here,  unlike  in 
the  case  of  shear-layer  modes,  the  intensity  of  the  oscillating  mode  is  independent  of 
background  disturbances.  The  streamwise  growth  rate  associated  with  the  oscillating 
mode  is  not  abnormally  large,  however.  The  frequency  of  the  oscillating  mode 
compares  well  with  predictions  based  on  a  spatio-temporal  theory,  but  not  with  those 
of  the  standard  spatial  theory. 

From  high-speed  films  it  is  found  that  the  overall  structure  of  the  oscillating  mode 
repeats  itself  with  extreme  regularity.  The  high  degree  of  repeatability  of  the  oscillating 
mode,  in  association  with  a  strong  pairing  process,  leads  to  abnormally  large  centreline 
velocity  fluctuation,  with  its  root-mean-square  value  being  about  30%  of  the  nozzle 
exit  velocity.  Energetic  and  highly  regular  pairing  is  found  also  to  lead  to  the  early  and 
abrupt  breakdown  of  the  potential  core.  The  regularity  often  extends  even  to  the  finer 
structure  immediately  downstream  of  the  breakdown.  An  attempt  is  made  to  explain 
these  special  features  both  in  terms  of  the  large-amplitude  vorticity  field,  and  in  terms 
of  the  theoretically  predicted  space  time  evolution  of  wave  packets. 


1.  Introduction 

Except  at  very  low  Reynolds  numbers,  a  laminar  shear  layer  with  constant  density 
is  unstable  to  arbitrarily  small  disturbances.  These  disturbances  evolve  by  selective 
amplification  as  they  convect  downstream  (Cohen  &  Wygnanski  1987).  The  migration 
of  vorticity  resulting  from  the  instability  forms  axisymmetric  or  helical  structures 
accompanied  by  large-amplitude  fluctuations.  Generally,  these  structures  interact  with 
one  another  and  contribute  significantly  to  entrainment  and  mixing  by  convective 
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!'i<ii  ri  i.  Power  spectral  density  in  y  helium  jet  ( .S  —  0.14).  showing  large  peaks  at  its 
suhharmonie.  and  several  harmonics.  The  hot  wire  was  located  along  the  jet  centreline  at  I  33 
diameters  dow  nstream  of  the  nozzle  exit.  I  he  l’rci|tiency  /„  of  the  oscillating  mode  does  not  depend 
on  the  spatial  position  in  that  region.  (  Taken  from  Sreenivasan  et  at.  IWt.i 


transport  <  Winant  &  Browand  1074).  and  to  the  generation  of  acoustic  noise  (Crighton 
1975:  Sarohia  &  M ussier  1977;  Kibens  19X0).  Very  often,  engineering  interest  arises 
when  the  nozzle  fluid  density  //,.  is  significantly  less  than  the  ambient  fluid  density.  t>  . 
The  kinematics  and  dynamics  in  the  near  field  of  variable  density  jets  is  the  focus  of 
the  present  study.  Before  discussing  its  specific  objectives,  it  is  helpful  to  review  the 
previous  literature  on  the  subject.  This  is  best  done  by  reviewing  the  experimental  and 
theoretical  studies  separately. 

1.1.  Experiment 

Several  early  studies  indicate  that  the  mean  tlovv  evolu  ion  in  the  near  field  can  be 
sensitively  dependent  upon  the  density  ratio.  S  =  />,.//>,.  Corrsin  &  Uberoi  (1949) 
measured  the  temperature  and  velocity  fields  for  various  nozzle  gas  temperatures  in 
heated  air  jets.  For  temperatures  corresponding  with  S  =  0.95  and  0.62.  they  found 
that  the  potential  core  terminated  at  x/D  x  5.  where  v  is  the  distance  along  the  axis 
and  D  is  the  nozzle  orifice  diameter.  When  S  was  decreased  to  0.49.  however,  the 
potential  core  terminated  at  x/D  3.5,  followed  by  a  relatively  sharp  drop  in  the 
centreline  temperature  and  velocity.  Beyond  the  potential  core,  the  jet  width  for  S  = 
0.49  was  approximately  twice  that  for  S  =  0.95.  Data  obtained  by  Landis  &  Shapiro 
( 1951 )  and  Sforza  &  Mons  ( 1978)  in  heated  air  jets,  and  by  Chriss  (1968).  Tombach 
( 1969)  and  by  Abramovich  et  al.  (1969)  in  heterogeneous  jets  confirm  that,  below  .V  ^ 
0.60.  mixing  processes  are  significantly  enhanced  in  the  near  field.  Abramovich  et  al. 
(1969)  measured  the  spread  rate  of  the  shear  layer  and  found  that  in  addition  to  the 
dependence  on  the  density  ratio,  spread  rates  were  also  sensitively  dependent  upon  the 
iet  Reynolds  number  and  upon  the  initial  shear-layer  thickness. 

There  is  evidence  to  suggest  that  the  mean  flow  variations  described  above  may  be 
linked  to  abnormally  intense  instability  modes  sustained  only  when  the  density  ratios 
are  small.  Smith  &  Johannesen  (1986)  obtained  power  spectra  using  microphones 
located  m  the  acoustic  far  field  of  helium/argon  jets  issuing  into  air.  In  pure  helium  jets 


Fiawi  2  (<t)  Instantaneous  digital  image  of  the  aerosol  distribution  in  the  axial  plane  of  a  seeded 
llow.  A  single  pulse  from  a  YAG  laser  illuminates  the  How.  .S'  --  0.29.  (  Taken  from  Kyle  I98X.)  (h) 
A  typical  realization  oT  (he  helium  jet  issuing  from  the  nozzle  into  ambient  air;  the  jet  spreads  \er> 
rapidly  beyond  a  reasonably  well-defined  streatmvisc  position.  The  figure  is  taken  from  Sreenixasan 
cl  ul.  ( 19X9).  where  details  are  given  of  how  this  picture  was  obtained ;  see  also  it 2  of  the  text  Re  - 
1400.  S'  -  0.14.  nozzle  diameter  --  6  mm. 


they  discovered  an  oscillatory  jet  instability  for  Mach  numbers  below  0.56.  Reminiscent 
of  the  mean  flow  results  of  Abramovich  ct  ul.  (1969).  this  instability  was  found  to 
depend  upon  the  precise  nozzle  contour.  The  density  dependence  of  jet  noise  for  higher 
velocities  is  reported  by  Hoch  cl  ul.  ( 1973)  and  by  Chan  &  Leong  ( 1973).  Kyle  ( 1986) 
and  Sreeni  vasan.  Raghu  &  Kyle  ( 1989)  examined  spectra  obtained  along  the  centreline 
of  low-speed  helium/air  jets,  and  found  that  the  jet  instability  of  the  type  obsetved  by 
Smith  &  Johannc.sen  is  sustained  in  jets  with  density  ratios  at  least  as  large  as  0  8  These 
authors  showed  that  the  oscillatory  behaviour  can  dominate  the  flow  in  the  entire  near 
field  (see  figure  1 )  Monkevvitz  ct  ul.  ( 1990)  have  obtained  similar  results  in  heated  air 
jets. 

The  observation  that  the  spectral  density  function  is  dominated  by  a  discrete 
frequency  spike  of  large  magnitude  and  its  higher  harmonies,  quite  unlike  the  situation 
in  constant  density  jets,  suggests  the  occurrence  of  a  new  phenomenon.  The  time 
records  used  for  calculating  these  spectra  are  so  long  that  many  wave  crests  pass  by  the 
tixed  hot-wire  position  during  a  single  time  record.  In  order  to  produce  such  ‘spiky' 
spectra,  there  must  be  relatively  little  random  variation  from  one  passing  structure  to 
the  next,  or  the  modes  must  preserve  their  phase  over  many  cycles.  Stationary  images 
in  the  near  field  of  variable-density  jets  reveal  that  the  spiky  power  spectra  are  indeed 
accompanied  by  extremely  coherent  ring  vortices  formed  along  the  jet  column,  as 
exemplified  by  figure  2(a)  and  other  photographs  presented  by  Kyle  ( 1988).  Sreeni  vasan 
<>/  ul.  (1989)  and  Monkevvitz  cl  ul.  (1989).  We  shall  denote  this  henceforth  as  the 
'oscillating  mode’.  Radial  profiles  of  the  mean  temperature  (Monkevvitz  ct  ul.  1990) 
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suggest  that  these  structures  may  he  unusually  elVecltve  in  convective  transport.  The 
variable-density  jet  can  undergo  a  catastrophic  spread  near  the  end  of  the  potential 
core,  as  seen  in  figure  2(h)  (see  also  Sreenivasan  a  at.  1989;  Monkewitz  at  at.  1989). 
This  spread  is  associated  with  the  radial  ejection  of  fluid  in  the  form  of  side  jets' 
(Monkewitz  et  a!.  1989).  which  may  be  due  to  the  growth  of  azimuthal  instabilities  in 
the  highly  strained  region  between  the  vortex  rings  ( Liepmann  1991 ).  At  the  end  of  the 
potential  core,  the  ring  structures  vanish  (figure  2).  Power  spectra  generally  show  that 
the  subharmonic  reaches  its  maximum  towards  the  end  of  the  potential  core  (Kyle 
1986,  1988\  suggesting  the  occurrence  of  vortex  pairing  (Ho  &  Huang  1982). 

In  addition  to  variable-density  jets,  several  other  flows  such  as  wakes  (e.g. 
Strykowski  &  Sreenivasan  1990)  and  countercurrent  shear  layers  (e.g.  Strykowski  & 
Niccum  1992)  show  similarly  well-ordered  structure  under  certain  conditions.  Partly  in 
an  effort  to  explain  these  observations,  the  theoretical  framework  of  space  time 
instability  of  fluid  flows  has  seen  a  resurgence  in  the  last  few  years.  Since  we  wish  to 
remark  on  the  applicability  of  the  theory  to  our  experiments,  it  is  necessary  to  review 
the  theory  briefly. 

1.2.  Theory 

Even  though  the  organized  state  that  we  are  trying  to  explain  here  is  highly  nonlinear, 
the  theory -which  purports  to  explain  the  onset  of  this  state  is  linear.  Unlike 
temporal  instability  of  the  basic  state  where  one  considers  the  growth  in  time  of 
spatially  periodic  perturbations,  or  the  spatial  instability  where  one  considers  the 
steady-state  response  of  the  basic  state  to  spatially  localized  but  time-periodic 
disturbances,  it  has  been  found  more  useful  to  study  the  response  of  a  quiescent  system 
for  t  <  0.  perturbed  impulsively  by  a  disturbance  of  the  form  -  4(.v)4'(/).  The  resulting 
disturbance  takes  the  form  of  a  wave  packet,  as  discussed  by  Sturrock  ( 1958).  Briggs 
(1964).  Gaster  (1968«.  b )  and  Huerre  &  Monkewitz  (1985).  among  others.  The  wave 
packet  may  be  rewarded  as  the  Green  function  from  which  the  complete  response  to 
any  distributed  perturbation  can  be  deduced. 

Using  closely  related  formulations,  Briggs  (1964)  and  Gaster  (1968</,  h)  have  shown 
that  if  the  wave  packet  evolves  over  a  long  period  of  time,  then  the  motion  at  every 
point  (x.  t)  in  the  physical  plane  which  results  from  a  single  pulse  pcnurbation  can  be 
associated  with  a  specific  wavenumber  k.  Thus  for  quite  general  conditions  (Bers  1983). 
along  each  trajectory  x/t  ---  constant  the  fluid  motion  is  eventually  dominated  by  a 
single  complex  wavenumber  A*.  which  simultaneously  must  satisfy  the  two  relations: 

<wj(A*)/<'Ar  =  0.  (1.1) 

hov(k*)/Tk,  =  x/t.  ( 1.2i 

We  may  view  (1.1)  and  (1.2)  as  determining  k*(x/t)  for  any  given  x/t.  Note  that  the 
functions  fur,  <».  are  just  the  real  and  imaginary  components  of  the  dispersion  relation. 

Our  intention  is  not  to  review  the  mathematics  which  leads  to  ( 1 . 1 )  and  (1.2).  but  to 
point  out  its  physical  content  to  the  extent  needed  here.  Each  wave  excited  by  the 
impulse  perturbation  immediately  starts  to  grow  (or  attenuate)  and  propagate  at  a 
distinct,  finite  rate.  In  due  course,  the  energy  associated  with  certain  of  these  modes  will 
tend  to  prevail  in  particular  locations.  To  see  which  modes  prevail  and  at  what 
locations,  let  us  imagine  several  observers  of  this  system,  each  moving  at  a  different 
speed  x/t.  At  first,  each  observer  may  notice  many  of  the  excited  waves  that  are 
propagating  at  his  speed,  but  one  wave  will  eventually  dominate  his  field  of  view  ;  it  is 
the  wave  wmch  simultaneously  has  an  associated  real  group  velocity  that  matches  the 
observer's  speed  (equation  ( 1 .2)).  and  has  a  transient  growth  rate  n>i  that  is  a  maximum 
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in  the  sense  oft!  I)  ll  several  wavenumbers  A  *i  \  n  satislv  (I  h  and  (IJr  the  one 
which  has  the  largest  growth  rate  will  dominate  Because  tlieie  is  .mb.  ,t  single 
dominant  mode  at  Inrite  lime  lor  each  v  t.  the  dislui  bance  w  ill  appear  to  each  ohsei \er 
as  a  discrete  frequenev  wav  drain  having  wavelength  Aj(  v  /)  and  Irequcncy  A  '  i.  and 
will  appear  to  grow  in  space  and  time  as  -  e\p|  (k*(  \  Mtv  4  '•■,(/>  ”  i 1\. 

In  certain  systems.  tire  wave  packet  will  evolve  in  such  a  wav  that  the  observer 
moving  with  speed  v  i  -  n  will  see  t Ire  movie  with  wavenumber  A*HH  grow  with  time, 
i.e.  that  f»,(A  *(()))  -■  t)  \s  pomtcvi  out  earlier,  this  movie  alone  will  eventual!;,  dominate 
the  observer's  field  of  view.  I  his  behaviour  has  been  analvsed  explicitly  by  Cia'ter 
(1968u).  Briggs  (1964)  and  others  (see  the  review  by  lluvrrc  <.V  Monkewit/  I'tqot. 
and  is  called  ' absolute  instability  It  has  further  been  shown  that  w hen  <•<,( A  '(<l)>  •  0. 
discrete  frequenev  oscillations  will  evolve  near  the  origin  not  only  in  response  to 
an  impulse  perturbation,  but  also  in  response  to  nearly  any  spatially  localized 
perturbation  with  arbitrary  time  dependence  (see  Briggs  1964.  52.3.5).  Often  in  real 
Hows,  natural  excitation  resembles  a  continual  random  process  that  is  localized  in 
space  (e.g.  at  the  system  boundary  ).  It  is  therefore  compelling  to  conjecture  that  if  such 
a  flow  resembles  in  some  approximate  sense  an  absolutely  unstable  system,  it  may 
support  oscillatory  instabilities  of  the  sort  described  earlier. 

Gaster  &  Duvey  (1968)  have  examined  the  wavepackct  which  evolves  from  an 
impulse  perturbation  m  the  inviscid  wake  In  this  study.  A  *(  v  /)  and  ^«(  A  * )  were 
calculated  for  a  range  of  x/t.  The  particular  wake  profiles  which  they  examined  were 
not  absolutely  unstable.  Huerrc  &  Monkewit/  (1985)  have  evaluated  the  impulse 
response  in  a  plane  shear  layer  between  two  parallel  streams,  and  found  for  certain 
countercurrent  flows  that  /->,< A  *( (>) )  >  0.  i.e..  if  perturbed  locally,  this  idealized  system 
will  eventually  develop  oscillations  with  frequency  '•/,(£*(()))  at  lived  points  surrounding 
the  space  time  origin.  Koch  1 1985)  has  examined  the  wake  for  a  broad  range  of  profiles 
and  found  that  the  wake  can  also  be  absolutely  unstable  Pav  ithran  &  Redekopp  ( 1989) 
have  studied  the  impulse  response  for  variaoic-density  shear  layers,  while  Monkewit/ 
&  Sohn  (19X6.  1988)  have  examined  heated  round  jets.  These  latter  authors  found  for 
a  family  of  density  and  velocity  profiles  that  the  jet  is  absolutely  unstable  when  S  is 
sufficiently  small. 


1.3.  Motivation  for  the  present  work  and  the  organization  of  the  paper 

This  experimental  investigation  is  aimed  at  obtaining  a  clearer  physical  understanding 
of  unstable  modes  which  evolve  in  the  near  field  (or  the  transition  region)  of  variable- 
density  jets.  A  key  feature  of  these  experiments  is  that  they  are  well  controlled,  in  that 
they  explore  within  the  ranges  covered  the  dependence  of  the  flow  on  each  of  the 
important  governing  parameters  separately,  while  the  others  are  fixed.  We  have  chosen 
conditions  for  which  the  flow  is  incompressible  and  buoyancy  effects  arc  negligible. 

We  had  earlier  (Kyle  1986;  Sreenivasan  et  al.  1989)  performed  an  experimental 
investigation  of  the  .V-dependcnce  of  the  oscillatory  mode,  and  found  that  at  some  jet 
speeds,  the  measured  onset  value  of  S  roughly  corresponded  with  the  theoretical  value 
for  which  <■;,(£*(()))  first  became  positive  Monkewit/  et  ai.  ( 1990)  found  this  behaviour 
in  heated  jets  as  well,  and  further  observed  that  the  Strouhal  number  generally  falls 
within  the  range  predicted  by  their  theory.  Detailed  quantitative  comparisons  between 
theory  and  observation  for  corresponding  values  of  governing  parameter  values  have 
not  been  made,  primarily  because  of  the  lack  of  experimental  control  already  discussed. 
An  added  motivation  for  the  present  study  is  to  compare,  where  possible  lot  a 
relatively  wide  parameter  range,  the  predicted  parametric  dependence  of  the  onset  of 
the  instability  with  the  observed  behaviour. 
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(he  rest  (4  the  paper  is  organized  as  loilows  Alter  a  discussion  m  k2  <4  the 
experimental  conditions,  we  discuss  in  O'  the  parameters  governing  the  let  hehnviom 
near  the  onset  of ’oscillating  mode  and  argue  that  n  slioukl  depend  essentially  upon  the 
Reynolds  number,  the  density  ratio,  the  ratio  ol  nozzle  diameter  to  shear-layer 
thickness,  and  background  disturbances.  In  si 4.  we  discuss  measurement'  loeusing  on 
obtaining  a  time-resolved  kinematic  description  ol  the  oscillating  mode  and  its 
subsequent  breakdown  process;  using  a  set  ol  high-speed  lilnis  obtained  tor  a  ra:  o| 

parameter  values,  in  conjunction  with  laser  Doppler  veloennetry  measurements  we 
study  the  growth  and  interaction  ol  large-amplitude  waves  associated  with  both  the 
oscillating  mode  and  normal  shear  layer  modes.  In  k  A  we  concentrate  on  teatures  ol 
the  power  spectrum  obtained  in  the  region  close  to  the  nozzle  exit,  where  the 
disturbances  are  small;  we  study  the  intensity  and  the  Irequency  ol  each  kind  ol  mode 
as  functions  of  'he  four  governing  parameters  and  ol  the  streamwise  distance, 
bmphasis  is  placed  on  the  behaviour  ol  the  oscillating  mode  ne.u  cutset,  and  on 
distinguishing  the  different  kinds  of  instability  existing  in  variable-density  lets.  W  here 
appropriate,  comparisons  between  experiment  and  theory  are  noted.  As  an  aid  to 
understanding  the  origin  of  the  oscillating  mode,  we  discuss  in  kb  its  response  to 
acoustic  forcing  and  other  changes  in  the  environment.  The  paper  concludes  with  k~ 
in  w  hich  the  nature  of  the  oscillating  mode  is  discussed  in  terms  of  both  x  orticit \  and 
the  space  time  instability.  The  Appendix  summarizes  certain  issues  involved  in 
interpreting  the  hot-wire  measurements  made  in  the  inhomogeneous  How  considered 
here . 


2.  Experimental  arrangement 

The  jet  facility  consists  of  a  stainless  steel  settling  chamber  with  screens  ol  varying 
gradation  and  a  capability  for  attaching  matching  nozzles  of  different  exit  diameters. 
The  nozzle  and  the  settling  chamber  are  movable  while  the  optical  table  surrounding 
the  nozzle  is  fixed,  making  measurements  at  various  streamwise  locations  relatively 
easy.  Details  of  the  facility  are  described  in  Stein  ( I  969).  Measurements  were  made  with 
nozzles  of  9  3  and  13.3  mm.  both  machined  from  aluminium  stock.  The  9.3  mm  nozzle 
(contraction  ratio  190)  is  designed  using  a  cubic  equation  for  the  contour  (Hussain  & 
Ramjee  1976).  and  the  1.3.3mm  nozzle  (contraction  ratio  93)  is  designed  with  an 
ASME  series  low  //  contour  (Bean  1971). 

Various  densities  of  the  jet  fluid  were  obtained  by  mixing  helium  and  air  in  desired 
proportions.  Before  being  mixed,  the  helium  and  air  flow  rates  were  metered  separately 
using  calibrated  rotameters.  The  calibrations  were  carried  out  separately  by  correlating 
the  rotameter  readings  with  the  volumetric  discharge  (J  through  the  nozzle.  Q  was 
calculated  using 

Q  k[!/) (24) 

where  D  is  the  exit  diameter  of  the  nozzle.  ( j.  is  the  potential  core  velocity  in  the  exit 
plane  (found  to  be  uniform  except  for  the  boundary  layer).  4*  is  the  displacement 
thickness  associated  with  the  boundary  layer  at  the  nozzle  exit.  The  displacement 
thickness  was  measured  directly,  as  described  in  §3.  l\.  was  determined  by  measuring 
the  pressure  drop  through  the  nozzle  with  a  Bara  iron  pressure  transducer  (Type  370H- 
10).  and  using  Bernoulli's  formula. 

A  few  critical  velocity  measurements  were  made  using  a  1  Si  laser-Poppler 
vclocimclcr  (I.DV)  operated  in  forward-scatter  mode:  an  aerosol  of  water  droplets  was 
used  for  seeding  particle.  A  hot  wire  (5  pm  diameter.  0.0  mm  active  length),  operated 
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Fkh  rf  3,  Experimental  configuration  used  for  obtaining  high-speed  motion  pictures  of  the  aerosol 
field  in  a  gas  jet.  Framing  speed  was  7000  f.p.s.  Rt  and  ft.  are  two  coneaxe  cylindrical  reflectors  w  Inch 
create  an  intense  plane  of  illumination. 


on  constant  temperature  mode  on  the  DANTEC  55M01  anemometer.  was  also  used 
to  measure  velocity  in  air  jets.  In  heterogeneous  regions,  the  hot  wire  cannot,  in 
general,  be  used  for  determining  the  absolute  value  of  the  velocity,  because  the  signal 
£is  a  function  of  both  velocity  u{l )  and  local  helium  concentration  iff),  i.e.  li  =  £(»(/). 
c(t)).  One  can  always  measure  the  velocity  inside  the  potential  core  where  the  fluid  is 
homogeneous,  but  such  probe  intrusion  can  significantly  alter  the  stability  cha¬ 
racteristics  of  the  oscillating  mode  (see  Sreenivasan  el  at.  1989).  Unless  otherwise 
stated,  the  hot  wire  was  positioned  outside  the  shear  layer  at  r  =  \D  +  50.  (t  is  shown 
in  the  Appendix  that  normalized  power  spectra  obtained  at  this  location  are 
approximately  the  same  as  the  normalized  power  spectra  of  a  obtained  inside  the 
potential  core,  as  long  as  the  disturbances  are  small.  Spectral  density  curves  were 
obtained  on  a  single  channel  HP  Spectrum  Analyser  (model  3561  A).  This  device 
displays  a  discretized  spectrum  using  400  equal  subintervals,  or  "lines'  to  span  the 
desired  frequency  interval.  Several  standard  spectral  windows  can  be  used  with  the  fast 
Fourier  transform  (FFT)  algorithm.  We  have  always  chosen  to  use  a  "flat  top’  window. 
A  few  of  the  spectral  measurements  to  be  reported  are  concerned  with  the  'width'  of 
a  given  peak  in  the  power  spectrum.  For  such  measurements,  the  bandw  idth  and  centre 
frequency  of  the  spectrum  analyser  display  were  adjusted  so  that  the  spectral  peak 
spanned  a  number  of  lines  which  was  of  the  order  of  10J.  This  practice  was  found  to 
minimize  the  effects  of  discretization  on  this  type  of  measurement. 

High-speed  motion  pictures  were  made  of  the  jet  by  illuminating  the  flow  with  a 
narrow  sheet  of  continuous  laser  light  (7  W)  in  a  plane  coinciding  with  the  jet  axis. 
Motion  pictures  were  obtained  using  a  Woliensak  'Fastax  WF6'  16  mm  motion 
picture  camera  with  400  ASA  movie  film.  In  order  to  capture  several  frames  within  one 
oscillation  cycle,  it  was  necessary  to  film  at  high  repetition  rates,  and  the  film  speed 
ranged  up  to  7000  s  This  required  intense  illumination  for  adequate  quality  of  the 
images  The  laser  slice!  was  therefore  formed  using  a  'multipass cell"  as  shown  in  figure 
3.  The  laser  light  was  reflected  back  and  forth  between  the  two  concave  cylindrical 
reflectors  creating  an  intense  plane  of  illumination  without  ever  passing  through  a 
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diverging  lens.  This  method  is  discussed  by  Long  et  at.  (19X3).  The  nozzle  fluid  was 
made  visible  by  seeding  it  with  a  polydisperse  aerosol  of  water  droplets  doped  with  a 
fluorescent  dye  (sodium  fluorescene).  The  droplets  were  generated  using  a  TSI  aerosol 
generator  (model  9306).  Instantaneous  images  were  obtained  using  a  pulsed  Nd :  Y  AG 
laser.  In  this  case,  the  laser  sheet  was  formed  using  a  concave  cylindrical  lens  to  diverge 
the  beam  in  the  axial  plane,  and  a  convex  cylindrical  lens  to  converge  the  beam  in  the 
horizontal  plane.  The  laser  had  a  power  density  of  2x  i()T  J  s  1  per  pulse  and  a  pulse 
duration  of  about  10  ns.  The  images  were  captured  on  a  Photometries  ‘IC200'  CCD 
array  with  a  1300  x  1035  pixels,  and  were  processed  using  in-house  software  (Prasad  & 
Sreenivasan  1990). 


3.  Governing  parameters 

Separated  free  shear  flows  of  heterogeneous  composition  can  potentially  depend 
upon  a  large  number  of  parameters.  However,  by  limiting  the  range  of  flow  conditions 
and  by  exploiting  justifiable  approximations,  one  can  construct  a  set  of  experiments 
whose  outcome  depends  upon  just  a  few  non-dimensional  parameters.  We  chose  at  the 
outset  to  limit  this  study  to  Mach  numbers  M  —  U„/a,  less  than  0.3;  here  a,  is  the 
sound  speed  for  the  ambient  air.  For  most  of  the  experiments.  M  was  substantially 
smaller.  Under  these  conditions,  the  Mach  number  effects  are  uniformly  small. 
Experiments  of  Bradshaw  (1966)  in  isothermal  air  jets  suggest  the  same  (see  also  §5.1 ); 
straightforward  scaling  arguments  applied  to  the  full  energy  equation  (Kyle  1991) 
suggest  that  this  is  so  also  for  the  variable  density  case.  Thus,  the  non-dimensional 
governing  equations  depend  upon  the  Reynolds  number  Re  -  U,.  D/i\ ,.  the  Richardson 
number  Ri  =  (/>,.—/■>*) Dg/(>,  U*.  and  the  Schmidt  number  Sc  =  r,  /DfU.  ,lir.  Here. 
DHv  air  is  the  diffusivity  of  helium  into  air.  g  is  the  acceleration  due  to  gravity,  and  i\, 
is  the  kinematic  viscosity  of  the  nozzle  fluid,  which  itself  can  be  calculated  using  Wilke's 
( 1 950)  approximation  based  on  the  statistical  mechanical  theory  of  Chapman  Enskog. 
In  the  present  experiments,  the  Richardson  number  Ri  ^  1.4x10  ’.  and  so  the 
buoyancy  effects  are  negligible  (Kotsovinos  1975).  Further,  although  Schmidt  number 
effects  may  not  be  neglected  a  priori ,  Sc  is  not  an  experimental  variable  in  this  study 
because  DHe  ajt.  is  essentially  independent  of  the  relative  concentrations  of  helium  and 
air  (Geankoplis  1972). 

3.1.  Inflow  boundary  conditions 

The  flow  everywhere  upstream  of  the  nozzle  exit  is  slow  and  has  uniform  composition. 
Reynolds  number  similarity  therefore  assures  us  that  the  non-dimensional  velocity 
distribution  at  the  nozzle  exit  is  a  unique  function  of  the  Reynolds  number.  Figures 
4 (a,  h)  show  radial  profiles  of  the  streamwise  velocity  V(r)/Uv  obtained  in  the 
separated  boundary  layer  for  the  13.3  and  the  9.3  mm  nozzles,  respectively,  over  the  Rc 
range  covered  in  the  experiments.  The  profiles  are  similar  and  agree  closely  with  the 
Blasius  profile.  For  all  profiles,  S*/b  was  within  ±4%  of  the  Blasius  value.  Here,  0  is 
the  momentum  thickness  defined  by 


(U(r)/Uv)[]-U(r)/Uv]dr. 


and  a'*  is  the  equivalent  plane  displacement  thickness  defined  by 


a* 


[1  —  U(r)/U,]dr. 


(3.1) 


(3.2) 
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Figurk  4.  Variation  of  streamwise  velocity  in  the  nozzle  plane  (v  —  0.06  mm)  at  various  Reynolds 
numbers:  O.  Re  =  16000;  □,  Re  =  7000;  x  .  Re  =  5000;  *.  Re  =  3500;  +.  Re  =  2700;  .  /?<■  = 

2100.  (nr)  U(r)/Ut.  for  ASME  nozzle.  D—  13.3  mm;  ( h )  l,’{r)/i\,  for  cubic  equation  nozzle.  />  = 

9.3  mm.  (c)  Fluctuation  amplitude  u'/Uv  for  ASME  nozzle;  I)  =  13.3  mm.  (d)  Fluctuation  amplitude 
u'/U,.  for  cubic  equation  nozzle;  D  =  9.3  mm.  Measurements  are  obtained  using  a  standard  hot  wire 
in  air  jets. 

The  upper  limits  of  integration  rn ,  signify  the  location  where  U(r)/l\,  —  0.1.  The 
longitudinal  velocity  fluctuation  intensity  profiles  u'(r)/Ue  (figures  4 (<•,</))  exhibit 
peaks  at  (| D  —  r)/()  =  3,  this  being  typical  of  jets  with  initially  laminar  shear  layers 
(Hussain  &  Zedan  1978).  Centreline  values  i/'( 0)/t/,,  at  various  speeds  L  j.  were  in 
the  range  0.08  0.19%  in  the  9.3  mm  cubic  equation  nozzle,  and  0.17  0.23%  for  the 

13.3  mm  ASME  nozzle.  Both  the  mean  and  fluctuating  velocity  data  indicate  that  the 
boundary  layers  separating  at  the  lip  of  the  nozzle  are  laminar  over  the  Re  range 
considered.  Because  the  separated  boundary  layers  are  laminar,  and  their  shapes  are 
essentially  invariant,  the  parameter  0  fully  characterizes  the  normalized  velocity 
profiles,  so  that  an  appropriate  non-dimensional  ratio  is  D/d.  Figure  5  shows  D/d 
values  calculated  from  the  above  profiles.  For  all  the  data  presented  here.  D/d  will  be 
reported  from  which  the  corresponding  value  of  Re  can  be  found  from  figure  5. 

The  jet  is  isothermal,  as  confirmed  by  direct  measurements.  The  density  ratio  S  = 
pr/p ,  is  therefore  just  equal  to  the  ratio  of  the  molecular  weight  of  the  helium/air 
mixture  to  that  of  air  since,  for  low  Mach  numbers,  the  effect  of  pressure  on  the  density 
is  negligible.  With  these  assumptions,  it  is  easy  to  show  (Kyle  1991 )  that  all  other  (non- 
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Figuri  5.  Momentum  thickness  0/D  as  a  function  of  the  Reynolds  number  Re.  Q.  13  ?  mm 
ASME  ft oz.zle:  O-  9-3  mm  cubic  equation  nozzle 

dimensionalized)  thermophysical  properties  of  the  nozzle  fluid,  such  as  the  kinematic 
viscosity  ratio  are  uniquely  specified  by  S. 

3.2.  Overall  parameter  dependence 

The  analysis  so  far  suggests  that  all  non-dimensional  characteristics  of  the  flow  should 
depend  primarily  upon  the  global  parameters  D/d.  S  and  Re.  A  major  goal  of  this 
paper  is  to  evaluate  the  dependence  of  the  oscillating  mode  of  these  three  global 
parameters.  As  we  are  interested  in  the  growth  of  small  disturbances  in  the  near  field, 
it  is  also  important  to  take  account  of  the  effects  of  the  time-dependent  disturbances 
imposed  on  the  flow.  These  are  due  to  acoustic  noise  and  to  vortical  disturbances 
originating  upstream  of  the  nozzle.  Such  background  disturbances  may.  under  certain 
circumstances,  strongly  affect  the  flow,  though  they  may  be  weak  (Cohen  &  Wygnanski 
1987;  §5  of  present  paper). 

Finally,  it  is  of  interest  to  note  that  in  the  limit  Re  -  x  and  M  -0.  the  system  of 
equations  and  boundary  conditions  governing  the  heterogeneous,  isothermal  flow  is 
analogous  to  those  governing  the  heated  air  jet  (Brown  &  Roshko  1974:  Kyle  1991). 
Whenever  possible  we  shall  compare  our  measurements  with  results  obtained  by 
Monkewitz  el  a/.  (1989,  1990)  and  Raghu  &  Monkewitz  (1991)  in  heated  air-jet 
experiments  for  similar  values  of  D/0  and  S. 

4.  Evolution  of  highly  organized  structures  in  the  near  field 

Before  proceeding  further,  it  is  helpful  to  have  an  overall  kincmatical  description  of 
the  transition  from  organized  to  disorganized  states  in  variable-density  jets.  At  the 
instant  shown  in  figure  2(a).  an  axisymmetric  wave  appears  close  to  the  nozzle  and  is 
followed  by  two  vortex  rings  in  the  streamwise  direction;  beyond  this,  the  nozzle  fluid 
necks  down  and  finally  'pinches  off',  leading  to  a  disorderly  structure  downstream.  A 
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0.29 
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32 
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75 
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Tabu  1 

.  A/0  and  L/0  as  measured  directly  from  the  films 

:  D  -  \ 

1 3.3  mm. 

Values  represent  an 

average  of  100  successive  eddies,  lets  dominated  by  the  oscillating  mode  arc  indicated  in  bold  print 


region  of  fluid  having  a  relatively  low  concentration  of  nozzle  gas  surrounds  the  ) 

disorderly  structure  and  extends  beyond  the  range  of  the  image.  In  order  to  study  the 

history  and  interaction  of  these  structural  elements,  high-speed  motion  pictures  have  > 

been  made  of  the  transition  region  of  the  jet  issuing  from  the  13.3  mm  nozzle  for  nine  i 

different  combinations  of  the  governing  parameter  S,  D/t)  and  Re.  The  density  ratios  ■, 

range  from  unity  (air  into  air)  down  to  S  —  0.14  (pure  helium  into  air)  while  D/t)  lies 
roughly  between  50  and  91  (see  table  1).  LDV  measurements  are  correlated  with  the 

films  for  the  same  flow  conditions.  The  results  of  this  effort  are  summarized  below.  * 

4. 1 .  General  observations 

In  certain  visual  respects,  the  evolution  of  large-amplitude  structures  was  found  to  be 

similar  for  all  nine  flow  conditions  considered.  Instability  waves  which  form  in  the  < 

initially  laminar  shear  layer  are  always  axisymmetric;  between  1  and  2.5  diameters 

downstream  of  the  nozzle,  these  waves  fold  or  "break"  and  quickly  roll  up  to  form  a 

train  of  vortex  rings,  the  centres  of  which  continually  move  downstream.  These  are  the 

vortex  ring  structures  seen  in  figure  2(a).  When  the  oscillating  mode  is  present,  each 

vortex  ring  is  involved  in  exactly  one  pairing  inside  the  viewing  region  (0  ^  x/D  5$ 

4.75).  whereas  in  jets  not  supporting  the  oscillating  mode,  the  interaction  is  less 
repeatable  and  a  vortex  ring  occasionally  passes  from  view  without  coalescence. 

The  pairing  process  itself  can  be  seen  in  the  sequence  of  ten  successive  movie  frames 
of  figure  6.  These  frames  advance  from  left  to  right,  with  the  left-most  frame  in  the 
bottom  row  immediately  following  the  last  frame  in  the  top  row.  By  comparing  frames 
1  and  6,  2  and  7,  and  so  forth,  it  is  seen  that  the  overall  structure  in  the  transition  region 
repeats  itself  to  a  remarkable  degree;  even  the  fine  structure  downstream  of  the 
potential  core  termination  repeats  for  many  cycles.  It  is  evident  that  the  processes 
involved  in  the  growth,  interaction  and  breakdown  of  vortices  all  evolve  temporally  in 
a  periodic  fashion. 

While  these  pictures  are  quite  instructive,  they  do  not  tell  the  whole  story;  flow 
visualization  experiments  in  constant  density  jets  also  show  some  repeatability  of 
structure  interactions  at  low  Reynolds  numbers  below  about  104  (Becker  &  Massaro 
1968;  Browand  &  Laufer  1975;  Davies  &  Baxter  1977).  However,  for  nozzle  Reynolds 
numbers  of  the  order  105  and  above,  there  is  much  less  repeatability  in  constant-density 
jets;  further,  the  initial  instability  of  the  laminar  shear  layer  could  be  both  axisymmetric 
and  helical  (Browand  &  Laufer  1975;  Drubka  &  Nagib  1981)  quite  unlike  the 
situation  here. 

The  differences  between  the  two  classes  of  jets  are  quantified  in  the  remaining 
subsections. 


4.2.  The  disturbance  wavelength  and  wave-breaking  length 
Two  lengthscales  which  characterize  the  disturbances  in  a  shear  layer  are;  (i)  the  wave¬ 
breaking  length,  L.  and  (ii)  the  disturbance  wavelength.  4.  These  are  defined  in  figure 
7.  Their  values,  as  deduced  from  the  high-speed  motion  pictures,  are  summarized  in 


Figure  6.  Successive  movie  frames  showing  the  aerosol  distribution  in  an  axial  plane  of  a  seeded  flow.  S  =  0.29;  D/ll  =  50;  Re  —  3025; 
D  =  13.3  mm.  The  sequence  proceeds  from  left  to  right  and  from  top  to  bottom.  Framing  rate  =  §/„.  The  horizontal  bands  in  the  photographs 
are  artifacts  of  using  the  multipass  cell  (see  figure  3). 
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Figure  7.  Typical  pattern  of  disturbance  evolution  (after  Becker  &  Massaro  1968). 

table  1 ,  where  the  data  have  been  non-dimensionalized  using  0.  Each  entry  represents 
the  average  data  on  100  successive  vortices. 

For  flows  which  do  not  support  the  oscillating  mode,  A/0  values  shown  in  table  1 
vary  about  a  mean  of  35.7  with  a  standard  deviation  of  about  2.3.  This  variation 
appears  random,  hence  A/0  appears  to  be  essentially  independent  of  D/0  and  5  in  the 
ranges  considered.  We  note  that  the  inviscid  linear  theory  of  Michalke  (1971)  indicates 
that,  at  D/0  v  69,  the  most  spatially  amplified  mode  for  S  =  1.0  has  a  wavelength 
A/6  *  33.  Furthermore,  the  analysis  indicates  that  A/0  is  only  weakly  dependent  upon 
D/6  and  S;  for  example  the  most  spatially  amplified  mode  at  D/0  as  69  for  S  =  0.50 
has  wavelength  A/0  «  34.  These  theoretical  results  are  in  good  quantitative  agreement 
with  the  data  shown  in  table  1  for  jets  which  do  not  sustain  the  oscillating  mode. 

For  flows  which  do  support  the  oscillating  mode,  A/0  differs  significantly  from  the 
values  just  discussed.  At  D/0  =  69.  S  is  reduced  from  0.48  to  0.29,  which  encompasses 
the  boundary  for  the  onset  of  the  oscillating  mode,  and  A/0  increases  from  35  to  44. 
This  suggests  that  the  characteristics  of  the  oscillating  mode  differ  from  other  shear- 
layer  modes.  Measurements  of  the  passage  frequency  presented  in  §5  will  show  that  this 
is  indeed  the  case. 

The  L/6  values  presented  in  table  1  show  a  marked  distinction  between  the  range 
of  wave-breaking  lengths  in  flows  which  sustain  the  oscillating  mode  and  those  which 
do  not.  For  D/0  =  50  and  60,  reductions  in  S  also  correspond  with  reductions  in  L/6 
(by  57%  and  48%,  respectively).  The  implication  is  that  the  instability  waves  break 
much  closer  to  the  nozzle  in  the  presence  of  the  oscillating  mode.  Physical  intuition 
suggests  that  the  wave-breaking  length  depends  upon  both  the  streamwise  amplification 
rate  of  the  instability  waves  as  well  as  their  ‘initial’  intensity.  In  §5  it  will  be  shown  that 
the  intensity  of  the  oscillating  mode  is  indeed  closely  associated  with  its  large  initial 
value. 
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Figure  8.  Histograms  of  wavelength  A/D  and  wave-breaking  length  L/D  for  various  S:  D/ti  =  69; 
Re  =  5500;  D  —  13.3  mm;  sample  size  =  100.  From  top  left  to  bottom  right;  S  —  0.66,  rr  =  0.08  (rr 
here  is  the  standard  deviation  for  the  histogram);  5  =  0.48,  rr  —  0.09;  S  =  0.29,  rr  —  0.04;  S  —  0.66. 
<r=0.14;  5  =  0.48,  (7  =  0.16:  5  =  0.29,  (7  =  0.05. 


4.3.  Orderliness  in  kinematics 

Figure  8  shows  histograms  of  A/D  and  L/D  for  three  different  values  of  S,  for  D/0  ^ 
69.  The  standard  deviation,  a ,  characterizing  the  uncertainty  in  the  measurement  of 
any  single  passing  structure  was  estimated  to  be  ±^D.  It  is  seen  that  rr  for  the 
wavelength  data  is  roughly  twice  for  S  =  0.66  and  0.48  than  that  for  S'  =  0.29.  Recall 
from  table  1  that  at  D/6  *  69,  the  oscillating  mode  dominates  the  transition  region  for 
S  =  0.29,  and  is  absent  at  S  =  0.66  and  0.48.  The  standard  deviation  for  the  wave¬ 
breaking  length  is  roughly  three  times  larger  for  S  —  0.66  and  0.48  than  that  for  S  = 
0.29.  These  results  indicate  that,  for  flows  dominated  by  the  oscillating  mode,  unusual 
regularity  exists  in  the  wavelength  and  in  the  location  where  the  wave-braking  process 
begins. 

In  figure  9,  the  locations  of  twelve  successive  eddies  are  followed  from  frame  to 
frame  for  S  =  0.66  and  S  =  0.29,  with  D/0  %  69.  The  trajectory  of  each  individual 
eddy  is  plotted  from  the  wave-breaking  point  to  the  point  where  coalescence  begins 
(see  figure  6).  Once  individual  vortices  begin  to  wrap  around  their  neighbours,  their 
location  can  no  longer  be  tracked.  Focusing  first  on  S  =  0.66,  it  is  seen  from  the  first 
point  of  each  trajectory  that  wave  breaking  occurs  randomly  in  time  (location  along 
the  abscissa)  as  well  as  in  space  (location  along  the  ordinate).  This  result  is  expected 
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Movie  tVames 


Movie  frames 

FKiliRi;  9.  Eddy  trajectories:  D/t)  =  69;  Re  =  5500;  I)  =  1 .3.3  mm.  Each  trajectory  begins  where  the 
wave  "  breaks'  and  ends  at  the  start  of  the  coalesc'd  ee  process.  Straight  lines  are  least-squares  fits,  [a) 
S  -  0.66.  (6)  .S'  =  0.29. 

from  the  histogram  of  L  shown  in  figure  8  for  the  same  conditions.  From  the  final  point 
of  each  trajectory  it  is  seen  that  the  coalescence  process  also  exhibits  significant  random 
variation.  For  S  =  0.29,  however,  both  the  roll-up  and  the  coalescence  are  quite 
regular  (see  also  figure  8).  Such  regularity  is  exhibited  in  all  of  the  jets  which  support 
the  oscillating  mode.  From  the  combined  facts  that  pairing  induces  complete  pinching 
off  of  potential  core  fluid,  and  that  this  process  occurs  with  extreme  regularity,  one 
might  expect  that  an  axial  profile  of  the  time-averaged  nozzle  fluid  concentration 
should  exhibit  a  very  sharp  drop  in  the  region  where  (bis  process  occurs.  This  is  in  fact 
the  case  (see  also  Chriss  1968;  Tombach  1969). 

Notice  that  the  framing  rate  in  the  films  was  adjusted  to  he  such  that  the  duration 
for  ten  frames  coincided  with  the  subharmonic.  This  orderliness  seen  in  figure  6  shows 
that  the  entire  flow  field  approximately  repeats  itself  after  every  subharmonic  period. 


I  I  M  2  v* 
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h'ltii  Ri  10.  Power  spectra  of  u  at  various  axial  locations:  S  =  0.29;  />///  =  50;  Rc  =  5025;  /)  =. 
13.5  mm.  The  corresponding  x/D  arc  marked  on  the  curves.  Spectra  are  calculated  from  time  traces 
obtained  using  LDV. 


4.4.  Centreline  velocity 

The  longitudinal  component  of  velocity  has  been  measured  at  various  streannvise 
locations  along  the  centreline  using  (he  LDV  set-up  mentioned  in  $2.  Possible  effects 
of  seeding  on  the  flow  were  investigated  by  monitoring  the  change  in  the  frequency  of 
the  oscillating  mode.  The  frequency  as  measured  by  a  hot  wire  outside  the  shear  layer 
never  differed  from  the  seeded  case  by  more  than  2%.  A  data  rate  of  10000  s  '  allowed 
us  to  construct  a  reasonable  time  trace,  from  which  the  spectra  shown  in  figure  10  were 
calculated  for  .S’ =  0.29  and  D/D-  50  In  figure  II.  the  spectral  intensities  of  the 
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x/D 

Figure  1 1.  Spectral  amplitudes  (6i  the  fundamental  and  the  subharmonic  seen  in  figure  10: 

O-  fundamental;  □.  subharmonic. 

primary  and  subhtirmonic  modes  are  shown  explicitly  as  functions  of  x/D.  Notice  that 
the  ordinate  in  figure  10  is  compressed  and  the  abscissa  is  expanded  (say.  relative  to 
figure  1 )  so  that  the  spikiness  of  the  spectral  peaks  is  not  apparent  at  a  cursory  glance. 

It  is  useful  to  correlate  the  spectral  intensities  in  figures  10  and  1 1  with  the  motion 
of  the  flow  structures.  At  x/D  =  0.6.  the  power  spectrum  exhibits  a  single  peak.  The 
frequency  of  this  mode  corresponds  with  the  average  local  passage  frequency  of  the 
initial  axisymmetric  wave.  This  is  the  fundamental  mode,  which  continues  to  grow  in 
intensity  along  the  centreline  even  as  the  waves  seen  in  the  films  break  and  roll  up.  The 
subharmonic  appears  in  the  spectra  at  around  x/D  —  1.  and  attains  its  maximum 
intensity  along  the  centreline  at  x/D  —  3.  It  is  seen  in  the  films  that  this  maximum 
occurs  at  a  location  where  fluid  along  the  centreline  is  rapidly  accelerated  during  the 
pinching  off  process,  which,  of  course,  occurs  once  every  subharmonic  period.  Beyond 
x/D  =  3  all  discrete  inodes  in  figures  10  and  1 1  continually  decay  in  the  streamwise 
direction  until,  at  x/D  ~  7,  the  subharmonic  is  barely  discernable  above  the 
background.  Correspondingly,  the  remnant  of  the  pairing  process  appears  in  the  film 
to  become  more  diffuse  as  it  travels  downstream  of  x/D  =  3. 

Figure  12  shows  the  centreline  profiles  of  the  root-mean-square  (r.m.s.)  velocity 
u'/Lj.  and  the  mean  velocity  V /  for  S—  1.0  and  0.29.  and  for  D/ti  50.  For 
reference,  note  that  these  arc  the  same  conditions  as  for  figures  6.  10  and  11.  The 
centreline  turbulence  intensity  for  S  =  0.29  attains  much  larger  amplitude  compared 
with  the  air  jet  (5  =  1).  To  our  knowledge,  such  large  intensities  of  the  order  0.3  do  not 
">ccur  in  an  unforced  constant-density  jet.  The  peak  u'/Uv  value  occurs  at  x/D  -  3.0. 
which  also  corresponds  with  the  peak  in  the  subharmonic  mode  (figure  1 1).  Thus,  the 
abnormally  large  velocity  fluctuation  found  in  this  jet  is  primarily  associated  w  ith  the 
vortex  pairing  process,  i  f  V„  is  seen  to  fall  oil  rapidly  starting  from  a  diameter  or  so 
from  the  nozzle.  On  the  other  hand,  we  have  already  mentioned  in  5j4.3  that  the  density 
of  the  gas  along  the  centreline  is  approximately  uniform  until  the  pinching  off  starts  at 
about  x/D  —  2.5.  In  the  region  l.0<.v/D<2.5  slow,  dense  ambient  air  is  being 
engulfed  by  the  roll-up  process,  but  this  mixing  does  not  penetrate  to  the  centreline. 
The  attenuation  of  V/V,.  therefore  implies  the  existence  of  an  adverse  streamwise 
pressure  gradient. 

In  summary,  the  breakdow  n  of  the  highly  organized  motion  near  the  nozzle  occurs 
in  a  qualitatively  similar  manner  in  all  jets  with  50  sj  D/0  ^  100  and  0.6  ^  S  s:  0.14. 
and  is  characterized  by  the  growth  of  axisymmetric  instability  waves  in  the  initially 
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Figuri;  12.  Streamwise  component  of  the  root-incan-squarc  ai  d  the  mean  velocity  as  a  function  of 

x/D:  O.  S  —  0.29:  □.  S=  1.0:  D/H=  50:  Re  =  3025 :  D=  13.3  mm. - ,  Corrsin  &  Uberoi 

(1049). 

laminar  shear  layer;  they  pair  and  pinch  off  the  potential  core.  The  unique  features  of 
the  oscillating  mode  are  the  following:  (i)  their  wavelengths  are  not  well  predicted  by 
linear  stability  theory  (§4.2);  (ii)  waves  break  very  close  to  the  nozzle  (§4.2);  (iii)  at  each 
point  in  the  near  field,  the  flow  is  highly  organized  around  the  subharmonic  period 
(§4.3);  and  (iv)  fluctuations  arc  intense  and  interact  strongly  with  the  mean  flow  field 
(§4.4) 

5.  The  dynamical  aspects  of  the  near  field 

In  this  section,  properties  of  the  oscillating  mode  are  studied  at  short  downstream 
distances  x/0  ^  75.  We  make  use  of  power  spectral  density  curves  obtained  from  a  hot 
wire  located  outside  the  shear  layer,  where  these  curves  are  representative  of  the  power 
spectrum  for  the  streamwise  velocity  inside  the  shear  layer  (see  Appendix).  These  ne..r- 
field  measurements  will  be  used  in  §§5.1.  5.2  and  5.3  for  the  purpose  of  studying  the 
evolution  of  the  oscillating  mode  with  respect  to  the  global  parameters  Re.  D/0  and  .S’. 
When  cast  in  dimensionless  form,  this  behaviour  may  be  meaningfully  compared  with 
the  theory  of  absolute  instability. 
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fua  Hi  13.  Power  spectra  ot'  u'  at  various  jet  velocities  (  \  :  \  /)  0.10  Spectra  are 

labelled  with  (  in  ms  ’.  The  arbitrary  relereitee  for  the  dB  scale  vane'  from  curve  to  curve. 

5.1.  Modal  variations  with  V  fixed 

We  emphasize  that  in  quantitative  studies  reported  in  this  paper,  one  parameter  at  a 
time  was  varied  keeping  others  lived.  The  first  section  is  meant  to  illustrate  the 
qualitative  changes  which  occur  in  the  power  spectrum  by  steadily  varying  the  jet 
velocity  while  holding  the  nozzle  gas  density  constant 

Figure  1 3  depicts  these  qualitative  changes  for  S  =  0.29  as  the  exit  v  elocity  is  varied 
from  3.0  to  9.0  m/s.  Note  that  the  spectral  distributions  corresponding  to  3.0  and 
5.0  m/s  are  similar  in  that  their  peaks  occur  at  identical  frequencies.  This  frequency 
invariance  was  also  exhibited  in  power  spectra  of  a  obtained  by  C  ohen  &  Wygnanski 
( 19X7)  inside  the  shear  layer  at  xj D  x  0  for  air  jets.  They  suggested  that  the  (relatively 
weak)  spectral  peaks  could  be  attributed  to  the  acoustic  resonance  properties  of  the 
upstream  plenum  chamber,  which  are  of  course  dependent  only  upon  the  cavity 
geometry  and  the  sound  speed  (see  also  Ctow  &  Champagne  1971).  The  same 
conclusion  had  been  reached  by  Hussain  &  Ramjee  (1976)  when  they  found  that 
neither  the  shape  of  the  nozzle  contour  nor  variations  in  D/ft  had  any  effect  on  the 
frequency  of  several  spectral  spikes  both  upstream  and  downstream  of  the  nozzle 
contraction.  These  observations  imply  that  the  distribution  of  spectral  peaks  obtained 
in  the  near  field  of  heterogeneous  jets  with  uniform  gas  composition  upstream  of  the 
nozzle  should  likewise  be  independent  of  the  jet  velocity,  as  indeed  evidenced  in  figure 
13.  We  also  note  that  the  frequency  of  fluctuations  that  are  excited  by  the  ubiquitous 
ambient  acoustic  signals  will  be  independent  of  the  velocity  and  density  . 

When  C,  is  increased  from  5.0  to  7.0  m/s.  however,  a  distinct  spectral  peak  emerges 
above  the  background  spectral  distribution  of  if.  This  is  the  oscillating  mode  which, 
near  to  its  onset,  does  not  appear  to  influence  the  intensity  of  modes  at  neighbouring 
frequencies.  Such  interaction  does  appear  only  when  the  oscillating  mode  has  grown 
more  intense  (as  in  the  uppermost  spectrum  in  figure  13  corresponding  to  i  — 
9.0  ni  .  s).  As  suggested  by  figure  13.  we  have  consistently  found  that  the  spectral 
intensity  of  the  oscillating  mode  increases  rapidly  and  continuously  with  I  \  near  1 1 
onset,  when  measuied  a*  anv  lived  location  (see  also  Monkcwitz  a  aJ.  1990).  The 
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Figirh  14.  Power  spectra  of  streamwise  velocity  fluctuations  u'  within  the  shear  layer.  S  -  0  14: 
v/ D  =  0.10.  Spectra  are  labelled  with  L  ,  in  ms  ’.  The  arbitrary  reference  for  the  dB  scale  varies  from 
curve  to  curve. 

frequency  of  the  oscillating  mode  increases  monotonicallv  with  Note  that  because 
the  oscillating  mode  alone  is  velocity  dependent  very  near  the  nozzle,  it  can  be  readily 
distinguished  from  the  background  disturbance  peaks  near  to  its  onset. 

Figure  14  shows  qualitative  changes  that  occur  in  the  spectral  distribution  of  shear 
layer  disturbances  as  the  velocity  is  increased  above  values  for  which  the  oscillating 
mode  has  already  achieved  a  relatively  large  intensity.  In  the  spectrum  corresponding 
to  V =  41.5  m/s  (at  the  bottom  of  figure  14),  the  tallest  peak  at  /  =  1733  Hz  is  the 
fundamental  of  the  oscillating  mode.  Its  subharmonic  and  several  higher  harmonics  are 
also  evident.  Between  (./,,  =  49.6  m/s  and  65.8  m/s  (second  and  fourth  spectra  from  the 
bottom)  the  velocity  evolves  from  a  stable  periodic  state  to  one  dominated  by 
disturbances  with  relatively  broadband  spectral  content.  Note  that  although  the  pure 
periodicity  is  lost,  these  relatively  broadband  disturbances  are  still  to  be  regarded  as  the 
oscillating  mode,  simply  because  the  frequency  and  intensity  change  smoothly  with  ( 
This  transition  has  been  discussed  by  Kyle  &  Sreenivasan  ( 1988.  1989)  and  Monkewitz 
et  al.  (1990).  In  this  same  velocity  range,  the  subharmonic  of  the  oscillating  mode 
diminishes  in  intensity  at  this  streamwise  location  and  finally  vanishes.  As  the  velocity 
increases  above  65.8  m/s  (fourth  spectrum  from  the  bottom),  the  oscillating  mode 
diminishes  in  amplitude  and  finally  disappears. 

At  about  57.7  m/s  (third  spectrum  from  the  bottom)  two  other  broadband  modes 
emerge  with  frequencies  which  differ  by  approximately  a  factor  of  two.  Their 
frequencies  are  apparently  not  related  to  that  of  the  oscillating  mode.  As  the  velocity 
increases  above  57.7  m/s.  these  modes  intensify.  It  may  be  thought  that  compressibility 
affects  the  disappearance  of  the  oscillating  mode  or  the  emergence  of  these  other 
broadband  modes,  since  M2  —  0. 126  for  the  uppermost  spectral  density.  However,  this 
is  almost  certainly  not  so  because  they  occur  at  A#4  =  0(10  2)  for  8’  =  0.29  and 
M~  =  0(10  :1)  for  S  =  0.48  (see  §5.2).  The  possibility  of  Mach-number  dependence 
will  therefore  not  be  further  mentioned. 
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Our  measurements  were  made  over  a  six-year  period  using  several  settling  chamber 
configurations,  four  nozzles  with  differing  contraction  contours  and  relate  to  the 
following  parameter  ranges;  0.14  <:  S  <  1.0:  40  <  D/d  125:  800  <  Re  <  15000.  On 
the  basis  of  data  obtained  from  films  (§4).  centreline  LDV  measurements  ($4)  and  near 
field  hot-wire  measurements  described  in  this  section  so  far.  we  conclude  that  the  only 
modes  that  are  sustained  in  our  facility,  which  are  unique  to  helium/air  jets,  are:  n)  the 
oscillating  mode,  along  with  its  subharmonic  and  higher  harmonics,  and  (ii)  the 
broadband  modes  evidenced  in  figure  14  at  higher  values  of  Did  and  Re. 

Monkewitz  et  at.  (1990)  have  reported  that  the  heated  air  jet  will  support  two  distinct 
oscillatory  instabilities  having  unrelated  frequencies,  which  they  call  'Mode  I'  and 
'  Mode  11'.  The  bulk  of  their  observations  pertain  to  ‘  Mode  11'.  We  note  that  kinematic 
features  of  their  ‘  Mode  II'  are  similar  to  those  of  the  oscillating  mode.  In  spite  of  the 
existing  uncertainty  about  the  parameters  governing  these  instabilities,  there  appears 
to  be  good  agreement  between  the  S  range  for  Mode  II'  (.Vs?  0.62)  and  for  the 
oscillating  mode  (5  =$  0.61).  as  well  as  between  the  Stu  =  I D/L  values  for  the  two 
modes  ( 5:  0.45).  The  oscillating  mode,  like  ‘Mode  II'.  was  determined  by  experiment 
to  be  axisymmetric.  From  these  comparisons,  and  others  to  be  developed  further  in 
$}5.2.  5.3.  it  is  reasonable  to  suppose  that  'Mode  II'  supported  in  heated  jets  and  the 
oscillating  mode  in  He/air  jets  are  essentially  the  same. 

Little  is  known  about  the  physical  nature  of  these  broadband  modes;  for  example, 
we  have  not  determined  their  azimuthal  dependence.  It  is  not  known  whether  they 
occur  in  heated  air  jets,  because  Monkewitz  el  al.  did  not  examine  L\.  beyond  where 
'Mode  IT  (the  oscillating  mode)  began  to  broaden,  whereas  it  is  precisely  at  these  l 
values  that  the  broadband  modes  first  appear  in  He/air  jets  (figure  14).  As  to  the 
possible  correspondence  of  these  modes  with  'Mode  L  of  Monkewitz  et  al..  it  should 
be  noted  that  ‘  Mode  I '  is  very '  spiky '  whereas  the  broadband  modes,  by  definition,  are 
not.  Furthermore,  at  any  given  .S'.  "Mode  I'  occurred  in  the  facility  of  Monkewitz  et 
al.  at  only  the  small  Ut.  values  (small  D/d  and  Re),  whereas  the  broadband  modes  were 
observed  in  our  facility  only  at  high  speeds  (large  D/d  and  Re).  We  tentatively  conclude 
that  the  'Mode  F  instability  does  not  occur  in  helium /air  jets. 

5.2.  D/d  dependence  of  the  frequency  of  the  oscillating  mode 
The  principle  of  dynamic  similarity  assures  us  that  when  fluctuation  frequencies  are 
expressed  in  non-dimensional  form,  their  values  must  depend  only  upon  D/d.  Re  and 
upon  background  disturbances,  if  .S'  is  fixed.  For  the  case  of  fixed  nozzle  gas 
composition  these  parameters  cannot  be  varied  separately  using  a  single  nozzle.  In 
order  to  examine  whether  the  Strouhal  number  .S7,,  =fD/Vv  is  a  function  of  either 
D/d  or  Re  alone,  data  obtained  in  both  the  ASMb  nozzle  and  in  the  cubic  equation 
nozzle  are  plotted  together  in  figure  15.  It  is  assumed  that  difference  in  the  exterior 
geometry  of  the  nozzle  lips  are  not  significant  enough  to  affect  the  results.  57, ,  values  for 
the  two  jets  are  quite  disparate  when  plotted  as  a  function  of  Re  (figure  15n).  where 
they  match  well  when  plotted  as  a  function  of  D/d  (figure  156).  Re  values  in  each  jet 
differ  by  approximately  a  factor  of  2.3  at  each  D/d  over  the  range  of  interest:  the 
maximum  initial  shear-layer  disturbance  levels  (»//(/,,)  for  fixed  D/d  values  differ  in  the 
two  nozzles  by  a  factor  of  about  2.2.  It  follows  that,  to  within  experimental  error,  non- 
dimensional  frequencies  for  the  oscillating  mode  as  well  as  for  other  modes  sustained 
at  large  D/d  are  dependent  upon  D/d  alone  over  the  parameter  range  of  interest,  and 
that  the  Reynolds  number  is  felt  only  indirectly  through  its  influence  on  D/d.  Subbarao 
(19X7)  studied  an  oscillating  mode  in  buoyant  helium  jets  (0.08  ^  Ri  fC  0.79)  having 
f  ully-developed  parabolic  velocity  profiles  at  the  nozzle  exit  ( D/d  =  15).  and  found  that 
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Fiouri;  15.  Non-dimensional  frequency  fl)/L„  plotted  separately  as  a  function  of  D/0  and  as  a 
function  of  Re:  0-0  =  9.3  mm:  □.  D  —  13.3  mm.  S  =  0.14:  x/D  =  0.10.  Solid  symbols  represent 
the  oscillating  mode:  open  symbols  represent  the  other  broadband  modes  prevalent  in  the  near-field 
measurements  (see  figure  14). 


Fiouri;  16.  Non-dimensional  frequency JD/l  ,.  as  a  function  of  D/D  for  various  values  of  S:  x/D  = 
0.10:  from  left  to  right  .S' =  0.14,  0.29  and  0.48  respectively  0.  oscillating  mode  from  present 
experiments:  O.  additional  prominent  disturbance  modes  in  present  facility:  @.  oscillating  mode 
(Mode  II)  in  the  heated  air  jet  (Monkewit?.  el  a!.  1990). 


.S7n  was  independent  of  Re  (1500  Re  <  12500).  However,  it  is  unclear  how  relevant 
his  results  are  to  the  current  ones,  as  the  parameter  ranges  for  the  two  studies  are 
different. 

In  figure  16.  the  non-dimensional  fluctuation  frequencies  shown  in  figure  15  for 
S  -  0.14  are  replotted  along  with  similar  data  obtained  for  -V  =  0.29  and  .S’  =  0.48. 
It  is  seen  that  as  .S'  changes,  the  functional  relation  between  the  near-field  stability 
characteristics  and  D/O  also  changes.  The  lower  onset  value  of  D/0  for  the  oscillating 
mode  is  about  40  for  S  —  0.2.9  and  about  60  for  S  —  0.48.  whereas  for  S  =  0.14  it 
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Figuri;  17.  Non-dimensional  frequency  fO/l  as  a  function  of  /)///:  x/D  =  0.10:  •  oscillating 
mode  from  present  experiments:  from  left  to  right  S  =  0.14. 0.29  and  0.48  respectively.  Q.  additional 
prominent  modes  in  present  facility;  ©.  oscillating  mode  in  the  heated  air  jet  (Monkewitz  ct  at.  1990); 
- -  theory  (Monkewitz  &  Sohn  1988). 


extends  below  40,  and  could  not  be  ascertained.  A  qualitative  feature  common  to  all 
three  density  ratios  is  that  the  oscillating  mode  vanishes  when  D/0  becomes  large 
enough.  At  least  one  of  these  new  modes  always  persists  beyond  the  D/0  bounds  for 
this  study.  Experimental  data  obtained  by  Monkewitz.  et  al.  (1990)  in  a  heated  air  jet 
with  S  —  0.47  are  plotted  along  with  helium/air  jet  data  for  .9  =  0.48.  D/0  ratios  were 
calculated  assuming  the  vorticity  thickness  S,„  ~  50  for  the  initial  velocity  profile  in  the 
heated  air  jet  of  Monkewitz  et  al.  (1990).  It  is  seen  that  the  comparison  is  very  good 
indeed,  even  though  at  each  D/0.  Re  differs  by  roughly  a  factor  of  2  and  the  two  flows 
are  not  strictly  analogous.  This  confirms  the  correspondence  between  ’Mode  11'  in 
heated  jets  and  the  oscillating  mode  in  helium/air  jets. 

The  broadband  modes  that  appear  when  D/0  reaches  large  values  one  or  two. 
depending  on  the  density  ratio  do  not  correspond  to  ‘Mode  F  of  heated  air  jets  in 
terms  of  the  frequency  range  (0.25  ^  Stu  ^  0.3)  as  well  as  density  ratio  (0.55  $  .V  s? 
0.69). 

The  frequency  data  shown  in  figure  16  have  been  non-dimensionalized  using  the 
momentum  thickness  as  the  lengthscaie,  and  the  resulting  St„—f0/l j,  values  are 
plotted  in  figure  17.  It  is  clear  that /scales  neither  on  D  nor  on  0  within  the  D/O  range 
covered.  In  order  to  find  out  whether  /  scales  'more  closely"  on  0  or  D.  we  have 
evaluated  A(.S7II)/(.S7l))tll,lx  for  each  of  the  curves  in  figure  16.  where  A(.S7,,)  is  the  range 
in  any  one  curve  and  (.SVI()m,lx  is  the  maximum  value  in  the  same  curve,  and  compared 
these  ratios  with  A(St„)/(Sr0),miK  calculated  from  figure  17.  At  each  S.  the  ratio  for  St„ 
is  approximately  half  as  large  as  for  57,,.  Further,  much  of  the  variation  in  57 ,,  occurs 
for  D/0  values  above  80  (figure  16),  whereas  St„  variations  arc  small  in  this  region 
(figure  17).  This  suggests  that  as  the  shear-layer  thickness  becomes  small  relative  to  the 
diameter,  the  shear-layer  thickness  tends  to  play  a  more  dominant  role  in  establishing 
the  frequency  of  the  oscillating  mode.  The  relative  invariance  of  St„  was  not  observed 
in  our  previous  study  (Kyle  &  Sreenivasan  1989)  because  of  the  limited  D/O  range 
considered. 

Another  important  feature  emerges  from  figure  17.  The  frequency  of  the  stationary 
mode.  rar(k*( 0)).  as  calculated  by  Monkewitz  &  Sohn  (1988).  is  compared  with  the 
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Figlki  18.  Near-field  basc-to-pcak  intensity  as  a  (unction  of  D,":  .  I)  -  1 5 . 3  mm: 

®.  D  -  9.3  mm.  S  =  0.14:  x/D  =  0.10.  Corresponding  Reynolds  numbe-s  are  shown  in  figure  15. 
(See  text  for  definitions.) 

experimental  data  for  5  —  0.48.  The  correspondence  seems  reasonably  good  for 
smaller  D/0,  but  degrades  as  D/0  is  increased.  The  Strouha!  number  of  the  broadband 
modes  also  agrees  with  the  theory  fer  a  range  of  D/0.  The  spatial  theory  predicts  that 
the  most  unstable  mode  for  5  =  0.5  should  have  0.01 25  C  St„  $  0.01 3  for  this  range  of 
D/0  (Michalke  1984).  which  is  clearly  a  poor  prediction. 

5.3.  The  intensity  of  the  oscillating  mode 

We  turn  our  attention  to  the  intensity  of  the  oscillating  mode.  Specifically,  we  are 
interested  in  choosing  a  non-dimensional  measure  of  intensity  which  will  reflect  the 
transition  that  the  oscillating  mode  undergoes  from  a  pure  periodic  disturbance  to  one 
with  a  relatively  broadband  spectral  content  (see  figure  14).  One  possible  measure  is  the 
ratio  of  the  maximum  intensity  on  the  spectral  peak,  located  al  /„.  so  the  intensity 
at  the  base  of  the  peak.  The  ‘base"  is  arbitrarily  chosen  to  be  wherevei  the  width 
of  the  spectral  peak,  measured  in  Hertz,  is  equal  to  -/„•  Thus,  letting  .S"„(  / )  denote  the 
power  spectral  density  of  1 T  at  frequency  f.  we  define  the  base-to-peak  intensity  as 
[SM/SM^.)]*.  This  measure  of  base-to-peak  intensity  is  generally  greatest 
whenever  velocity  fluctuations  associated  with  the  oscillating  mode  arc  simultaneously 
intense  relative  to  background  disturbance,  and  distributed  over  a  spectral  bandwidth 
smaller  than 

In  figure  18.  the  near-field  base-to-peak  intensities  for  jets  issuing  from  each  of  the 
nozzles  are  plotted  as  functions  of  D/0  for  S  —  0. 14.  The  curves  correspond  to  fourth- 
order  polynomials  fitted  to  the  data.  The  curve  for  the  ASM  E  nozzle  matches  quite  well 
the  curve  for  the  matched  cubic  nozzle.  When  the  data  are  plotted  as  a  function  of  Re. 
the  two  curves  are  relatively  far  apart  (this  result  is  not  shown).  Thus,  given  S.  the  near- 
field  base-to-peak  intensity  is  mainly  dependent  upon  D/0  alone.  We  had  formerly 
remarked  (Srcenivasan  et  al.  1989)  that  the  transition  to  an  oscillatory  state  occurs 
abruptly  with  respect  to  changes  in  the  governing  parameters.  This  is  illustrated  by  the 
rapid  rise  in  base-to-peak  intensity  as  D/0  increases  from  its  lower  limit  (see  also  figure 
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Figure  19.  Near-field  base-to-pcak  intensity  as  a  function  of  D/0  and  .S':  contours  are  labelled  with 
base-to-peak  intensity  values.  The  inner  contours  correspond  to  larger  and  larger  peak  intensities. 
Error  bars  denote  uncertainty  in  experimental  determination  of  onset.  ©.  onset  in  heated  air  jets 
(Monkewitz  et  al.  1990). 

13).  As  D/6  nears  its  upper  limit,  the  base-to-peak  intensity  varies  relatively  slowly. 
This  latter  aspect  is  reflected  in  Sreenivasan  et  al.  (1989)  by  the  large  error  bars  used 
by  them. 

In  summary,  we  have  shown  that  the  non-dimensional  frequency  of  the  oscillating 
mode  and  its  near-field  base-to-peak  intensity  are  definite  functions  of  S  and  D/6 
alone.  Both  these  quantities  vary  significantly  with  D/6  and  do  not  scale  with  either  /> 
or  6  alone.  An  upper  D/6  limit  for  the  oscillating  mode  has  been  found  for  all  S  0.14. 
This  limit  is  not  influenced  by  either  the  Mach  number  or  the  Reynolds  number. 
Indeed,  figure  1 5  shows  that  the  upper  D/6  bound  is  nearly  the  same  for  both  nozzles. 
A  lower  D/6  limit  also  exists  (except  perhaps  for  S  —  0.14).  Such  limits  are  useful  for 
assessing  the  accuracy  of  physical  theories  (see  §7)  which  predict  ‘critical  values'.  They 
are  also  important  in  the  practical  sense  that  they  help  establish  the  parameter  region 
relevant  to  the  appearance  of  the  oscillating  mode. 

5.4.  The  onset  of  the  oscillatory  state 

In  this  section  we  study  changes  in  the  flow  which  result  from  small  changes  in  ,5  near 
the  onset  value  S„  for  the  oscillating  mode,  for  various  fixed  values  of  D/6.  Physically, 
this  means  that  the  normalized  mean  velocity  profile  at  the  exit  is  held  fixed  as  the 
nozzle  fluid  density  is  changed  by  small  increments. 

Figure  19  shows  the  results  of  systematic  measurements  of  the  near-field  base-to- 
peak  intensity  for  the  range  of  ordered  pairs  ( S.D/6 )  relevant  to  this  study.  To 
construct  these  contours,  base-to-peak  data  were  obtained  for  S  =  0.14.  0.29.  0.4S  and 
0.60.  At  each  .S’,  curves  were  fitted  to  the  data  exactly  as  for  figure  18.  Data  obtained 
from  tiicse  curves  were  used  to  construct  the  contours  in  figure  19.  Onset  values  were 
further  explored  by  slowly  decreasing  S  at  fixed  D/6  values  of  83.  70.  60  and  50.  Error 
bars  are  shown  for  onset  values  only.  The  outermost  curve  marked  ‘O'  is  the  locus  of 
the  onset  value,  .S’„,  for  the  oscillating  mode  and  the  corresponding  D/6.  The  region 
inside  this  curve  (i.e.  below  and  to  the  left)  correspond  to  inlet  flow  conditions  which 
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support  the  oscillating  mode.  Everywhere  outside  this  curve,  the  oscillating  mode  is 
absent.  The  S0 ,  D/6  curve  is  useful  for  comparing  the  onset  of  the  oscillating  mode  with 
the  predicted  onset  -  or  ‘critical’ -- values  of  S  and  D/6  associated  with  absolute 
instability.  Over  a  range  of  D/6  values  which  includes  D/6  —  83,  critical  S  values  are 
predicted  to  lie  between  0.66  and  0.72  (Monkewitz  &  Sohn  1988).  Thus,  in  the  vicinity 
of  D/6  =  83  (S„  =  0.61),  one  might  say  that  there  is  reasonable  correspondence  with 
the  theory.  On  the  other  hand,  consideration  of  S„  over  a  wider  range  of  D/6  values 
leads  to  a  salient  discrepancy:  while  figure  19  shows  that  for  each  S  ^  0.61  there  is  an 
upper  bound  on  the  D/6  values  for  which  the  oscillating  mode  can  exist,  the  theory 
predicts  that  for  all  S'  <  0.62,  «j(/;*(0))  >  0  even  in  the  limit  of  D/6-  •*  x  (Monkewitz 
&  Sohn  1986,  1988). 

We  are  interested  in  the  transition  to  a  stable  periodic  state.  Note  that  the  contours 
in  figure  19  are  labelled  by  the  base-to-peak  intensity.  By  drawing  horizontal  cuts 
through  these  contours  at  different  constant  values  of  D/6,  one  immediately  sees  that 
the  transition  to  a  stable  periodic  state  brought  about  by  decreasing  S  can  occur  in 
various  ways,  depending  on  the  choice  of  D/6.  For  example,  at  D/6  =  50.  a  small 
decrease  in  S  below  So«0.31  is  accompanied  by  a  dramatic  rise  in  base-to-peak 
intensity,  whereas  at  D/6  =  83  base-to-peak  intensity  varies  slowly  with  S  below  S„  ^ 
0.61. 

In  order  to  quantify  the  effects  of  D/6  further,  the  square  root  of  the  power  spectral 
density  of  u'  at /0,  uf  ,  is  shown  in  figure  20  non-dimensionalized  using  t/t„  for  D/6  — 
50,  62.5,  and  83.  The  data  are  plotted  as  functions  of  the  normalized  density  difference 
AS  =  (Sn  —  S)/SQ.  The  spectra  were  measured  using  a  hot  wire  located  at  ,v  =  10ft.  and 
r/D  =  0.55  for  all  three  D/6  values.  Shadowgraph  images  for  each  flow  condition 
revealed  no  detectable  helium  present  at  these  locations. 

Figure  20  shows  that  the  functional  relation  between  u’f  JUv  and  AS  is  different  at 
different  values  of  D/6.  The  uncertainty  in  AS  results  primarily  from  the  uncertainty 
in  the  empirically  determined  value  of  S„.  For  small  AS,  the  data  in  figure  20  may  be 
fitted  to  power  laws  of  the  form : 

ufJ U(,  =  a(  AS)"  (AS  >  0),  j 
u/JU,  =  0  (AS<0),J 

where  the  free  parameters  a  and  n  are  determined  by  the  fitting  algorithm.  For  D/ft  = 
50,  a  =  101  ’  and  n  =  1.3,  while  for  D/6  =  62.5,  a  =  1 0* "  and  n  —  1.36.  For  D/ft  =  83. 
the  gradient  changes  snarply  at  AS  5;  0.15,  so  that  significant  error  occurs  if  the  power 
law  relation  (5.1)  is  fitted  to  the  data  spanning  the  entire  interval  0  ^  AS  $  0.76.  If  the 
model  (5.1 )  is  applied  only  to  points  lying  inside  the  interval  AS  =$  0.15.  however,  the 
power-law  assumption  seems  suitable,  with  a  =  10  1:1  and  //  =  0.55.  It  should  be 
mentioned  that  the  measurements  were  repeated  at  two  other  radial  locations  with  the 
results  that  the  coefficient  a  becomes  smaller  for  points  further  away  from  the  jet  (and 
is  thus  a  function  of  position),  but  the  exponent  n  was  essentially  independent  of  the 
probe  position. 

Monkewitz  el  al.  (1990)  measured  the  amplitude  of  near-held  pressure  fluctuations 
in  a  parameter  region  quite  near  to  the  latter  truncated  range  (onset  occurred  at 
D/ft  %  74,  S  —  0.62).  A  power-law  relation  provided  a  good  fit  when  they  assumed 
the  power-law  exponent  n  —  0.5  for  data  spanning  the  interval  AS  ^  0.26.  They 
further  noted  that  Landau's  weakly  nonlinear  stability  theory  (Landau  &  Lifshitz 
1959)  predicts  the  r.m.s.  amplitude  of  a  discrete  temporal  instability  should  lie  on  the 
parabola  u  l  /  AS  when  AS  is  small.  Figure  20  shows  that,  while  Landau's  theory  may 
apply  to  the  data  at  D/ft  =  83,  it  fails  for  D/ft  =  50  and  62.5. 
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Figure  20.  Square  root  of  the  power  spectral  density  of  u'  at  /„,  uj,  normalized  by  6',..  and  plotted 
as  a  function  of  AS  =  (S„  —  S)/S„  for  various  D/0:  The  power  spectrum  was  measured  with  respect 
to  an  arbitrary  reference  voltage  which  was  constant  for  all  of  the  measurements.  <>.  D/0  -  83: 
D/0  =  62.5;  O,  D/0  -  50.  S  =  0.14;  r/D  —  0.55.  The  uncertainty  in  S,  (see  figure  19)  is  the  primary 
source  of  uncertainty  in  AS. 

This  failure  is  not  unexpected.  Recall  that  Landau  &  Lifshitz  (1959)  derived  the 
equation  beginning  with  the  factorization 

=  A(t;  Sn-S)f(x).  (5.2) 

which  assumes  that  the  spatial  eigenfunction  f  is  independent  of  the  control  parameter 
S  near  the  critical  value  Slt.  The  data  of  figure  12  (see  §4.4)  suggest  that  Landau's 
assumption  concerning  f(x)  near  onset  does  not  apply,  say,  at  D/0  =  50.  The  large 
intensities  shown  in  figure  12  for  S’  =  0.29  correspond  with  AS  0.07.  which  may  be 
considered  small.  Even  so,  the  velocity  profiles  shown  in  figure  12  imply  that  the  mean 
flow  divergence  is  profoundly  altered  by  the  onset  of  the  oscillating  mode.  Recall 
further  from  table  1  (§4.2)  that  for  D/0  —  50  and  60,  small  reductions  in  S  correspond 
with  relatively  large  reductions  in  the  wave-breaking  length.  These  results  suggest  that 
the  mean  flow  in  the  transition  region  can  be  significantly  affected  by  the  onset  of  the 
oscillating  mode,  especially  at  D/0  of  the  order  50.  Figure  20  shows  that  near  onset, 
the  near-field  spectral  intensity  and  its  gradient  are  largest  at  just  these  D/0  values.  One 
may  therefore  conclude  that  for  D/0  values  of  the  order  50,  the  rapid  increase  in  near¬ 
field  disturbance  intensity  close  to  the  onset  can  lead  to  significant  variations  in  the 
mean  flow  structure.  If  the  spatial  distribution  of  fluctuation  intensity  is  related  to  the 
structure  of  the  mean  flow,  it  follow  that  (5.2)  should  not  hold  for  these  D/0  values. 

Although  there  is  considerable  variation  in  the  character  of  the  onset  of  the 
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oscillating  mode,  one  may  interpret  the  occurrence  of  a  stable  periodic  state  in 
homogeneous  round  jets  as  a  ‘bifurcation'  of  a  nonlinear  system  from  an  equilibrium 
state.  Although  there  are  obvious  problems  with  this  interpretation,  such  as  defining 
the  equilibrium  state,  the  bifurcation  is  a  supercritical  one  for  flow  conditions 
pertaining  to  this  study.  This  result  is  in  contrast  with  the  findings  of  Sreenivasan  et  at. 
(1989)  for  a  small  range  of  parameters;  they  used  a  4.0  mm  nozzle  and  found  that  the 
oscillating  mode  was  conditionally  stable  i.e.  that  the  transition  to  a  stable  periodic  state 
should  be  modelled  as  a  subcritical  bifurcation.  We  believe  that  this  difference  in 
behaviour  is  due  to  changes  in  the  mean  velocity  profiles  at  the  nozzle  exit ;  however, 
because  of  the  small  size  of  that  nozzle,  we  could  not  measure  the  velocity  profile. 

5.5.  Streannvise  evolution  of  the  oscillating  mode 

In  §4,  the  growth  and  interaction  of  large-amplitude  vortical  structures  were  examined. 
We  were  able  to  distinguish  kinematically  between  the  waves  associated  with  the 
oscillating  mode  and  waves  associated  with  shear-layer  modes  usually  studied  in 
homogeneous  jets.  In  this  section  we  aim  to  explore  this  distinction  further  by 
examining  their  streamwise  dependence  near  the  nozzle,  where  the  disturbances  are  still 
small.  We  first  establish  a  method  for  determining  the  most  spatially  amplified  mode 
between  any  two  streamwise  locations  in  the  air  jet,  and  then  extend  it  to  study- 
heterogeneous  jets.  For  reasons  that  will  be  made  clear  later,  this  method  is  most 
instructive  when  the  oscillating  mode  is  relatively  weak.  Therefore  most  of  the 
measurements  are  obtained  at  D/6  values  near  83,  for  which  it  has  been  shown  (see 
§4.2)  that  oscillation  intensities  are  small. 

Figure  21(«)  shows  power  spectra  for  u  measured  using  a  hot  wire  located  in  the 
shear  layer  of  an  air  jet  where  U/U(,  —  0.60  at  streamwise  locations  x/6  =  50.  100.  200 
and  300;  D/6  =  93  for  this  flow.  All  four  spectra  are  normalized  by  the  maximum 
spectral  amplitude  at  x/6  —  300.  It  is  evident  that  the  mode  at  670  Hz  plays  a  dominant 
role  in  the  evolution  of  the  shear  layer  for  these  conditions.  This  mode  is  prominent 
at  x/6  —  50  and  grows  continually  until  it  dominates  the  spectrum  at  x/6  =  200.  At 
x/6  =  300  the  subharmonic  has  become  dominant. 

We  would  like  to  characterize  the  growth  of  small  disturbances  in  a  region  of  the  jet 
where  nonlinear  effects  are  small.  With  this  intention  we  focus  on  the  region  x/6  <  100. 
Figure  21  ( b )  shows  the  ratio  of  spectral  amplitudes  at  x/6  =  100  to  those  at  x/6  =  50 
at  corresponding  frequencies/  This  function  simply  equals,  at  any  given  frequency,  the 
spatial  growth  rate  of  disturbances  integrated  between  x/6  =  50  and  x/6  =  100.  The 
mode  that  is  most  spatially  amplified  corresponds  to  the  maximum  value  of  this 
function. 

Comparison  of  the  total  amplification  function  (figure  21  b)  with  the  spectra  in  figure 
21  (a)  reveals  that  the  most  energetic  mode  at  x/6  -  200  (670  Hz)  is  not  the  most 
amplified  mode  in  the  near  field  ( a  860  Hz).  The  spectra  indicate  that  the  mode  at 
670  Hz  becomes  dominant  owing  to  the  combined  effects  of  its  near-field  prominence 
and  its  subsequent  amplification.  That  the  initial  velocity-independent  spectral 
distribution  of  u'  could  affect  the  frequency  of  the  most  energetic  mode  further 
downstream  in  the  shear  layer  was  shown  by  Cohen  &  Wygnanski  (1987);  see  also 
Becker  &  Massaro  (1967)  and  Gutmark  &  Ho  (1983).  The  total  amplification  curve 
(figure  216)  allows  one  to  distinguish  between  amplification,  which  varies  relatively 
smoothly  with  frequency,  from  the  influence  of  initial  field,  which  varies  randomly. 

Similar  information  can  be  extracted  from  the  total  amplification  function  in 
heterogeneous  jets,  as  long  as  the  disturbances  are  sufficiently  small.  In  figure  22.  the 
amplification  of  shear-layer  disturbances  between  x/6  =  25  and  75  is  shown  as  a 
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Figure  21.  (a)  Normalized  power  spectra  of  u'  obtained  in  the  shear  layer  of  an  air  jet  at  various 

streamwise  locations:  — .  x/O  -  50; . .  x/0  =100;  - - ,  x/O  =  200;  .  x/tt  =  300.  S  = 

1.0;  U/  Gj,  =  0.6;  D/0  =  93.  ( h )  Mode  amplification  between  50 0  and  1000  as  a  function  of  frequency 
using  spectra  shown  in  (a). 
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Figure  22.  Mode  amplification  between  x/O  =  25  and  x/O —  75  as  a  function  of  frequency  for 
various  values  of  .S'.  Curves  arc  labelled  with  S.  The  gain  at  each  .Vis  determined  to  within  an  arbitrary 
constant  factor.  Arrows  indicate  the  frequency  of  the  oscillating  mode. 
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Figl  ri  23.  Non-dimensional  frequency. fO/U,.  plotted  as  a  function  of  .9  for  (a)  D/0  =  83 :  (/>)  D/0  — 
73:  (<•)  D/0  =  62.5.  O.  oscillating  mode:  O-  most  amplified  mode  as  measured  between  x/0  =  25 
and  75.  Errors  bars  denote  measurement  uncertainty. 

function  of  frequency  for  S  values  ranging  from  1 .0  to  0.50.  for  D/0  —  8.3.  (For  reasons 
described  in  the  Appendix,  these  curves  arc  known  only  to  within  an  undetermined 
constant.)  Notice  that  the  curves  resemble  one  another  in  shape,  and  that  the  frequency 
of  the  most  amplified  mode  varies  smoothly  with  S.  This  continuity  suggests  that  the 
underlying  physical  processes  involved  in  the  streamwise  amplification  of  disturbances 
in  homogeneous  shear  layers  are  not  drastically  different  from  those  in  the 
inhomogeneous  shear  layer  for  this  range  of  density  ratios.  This  result  is  substantiated 
by  the  general  agreement  found  in  (54.2  between  the  wavelength  of  regular  shear-layer 
waves  and  the  theoretically  predicted  value.  To  our  knowledge,  this  result  has  not  been 
experimentally  confirmed  previously  (see  discussion  in  Michaike  1984).  Regarding 
the  oscillating  mode,  recall  from  figure  19  that  for  D/0  =  8.3  the  jet  supports 
oscillations  whenever  S  <  0.61 .  The  frequency  of  the  oscillating  mode  for  5  =  0.60  and 
0.50  is  indicated  in  figure  22  by  an  arrow.  We  note  that  the  frequency  of  oscillations 
is  distinct  from  that  of  the  most  spatially  amplified  mode.  Evidently,  for  D/O  =  8.3  the 
jet  supports  both  modes  simultaneously.  We  shall  shortly  see  how  this  observation  is 
consistent  with  the  observed  dominance  of  the  oscillating  mode. 

In  figure  23  the  non-dimensional  frequency  of  both  the  oscillating  mode  and  the 
most  spatially  amplified  mode  are  plotted  as  functions  of  5  for  D/O  —  83.  73  and  62.5. 
Notice  that  the  curves  for  the  two  modes  do  not  overlap  for  D/O  =  73  and  62.5.  For 
these  D/O,  the  oscillating  mode  is  very  intense  close  to  the  nozzle  even  for  S  near  S„. 
As  a  result,  the  harmonics  of  the  oscillating  mode  grow  rapidly  in  the  region  25  s; 
x/0  s$  75,  and  indeed  dominate  the  amplification  curve  for  S  <  Sn.  For  all  three  D/O. 
there  is  no  obvious  relationship  between  the  frequencies  of  the  most  spatially  amplified 
mode  and  of  the  oscillating  mode.  This  result  could  not  have  been  obtained  from  any 
scaling  argument,  simply  because  no  experimental  study  has  been  conducted  which 
shows  how  the  most  spatially  amplified  mode  varies  as  a  function  of  either  D/O  or  S. 
In  fact  it  is  clear  from  figure  23  that  direct  measurements  are  required  to  show  that  the 
frequencies  are  unrelated. 

We  note  that  fO/U,.  for  air  jets  shown  in  figure  23  are  close  to  those  measured  by 
Michaike  (1971),  but  are  some  30%  below  Drubka  &  Nagib's  (1981)  data  ( =  0.013). 
The  la'ter  authors  non-dimensionalized  the  frequency  using  0  measured  away  from  the 
nozzle.  Because  the  momentum  thickness  increases  with  v,  this  procedure  yields  a 
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Fiourh  24.  (a)  Mode  amplification  as  a  function  of  frequency:  .V  =  0.50;  D/0  =  S3.  Curves  are 
labelled  with  the  appropriate  streatmvise  interval.  The  gain  at  each  S  is  determined  to  within  an 
arbitrary  constant  factor.  (/>)  Power  spectra  of //' obtained  in  the  shear  layer:  .  x/tt  -  1; 

-v/W  =  75.  S  =  0.50:  D/0  —  <53.  The  arbitrary  reference  tor  the  dB  scale  dillers  between  the  wo  curves. 

higher  value  for  j'0/U,.:  Drubku  &  Nagib  indeed  noted  that  0  increased  by  40%  in  the 
first  80  momentum  thicknesses  downstream.  Regarding  the  theory,  the  maxim’. m 
spatial  growth  rate  is  predicted  to  occur  at  slightly  lower  t'0/V  values  as  S  is  decreased 
(Michalke  1984).  This  agrees  qualitatively  with  the  dan.  depicted  in  figure  23  It  should 
be  noted  that  the  absolute  value  ff*/Uv  corresponding  with  the  spat-ally  mod  amplified 
mode  was  found  to  be  some  30  40%  below  the  predicted  value,  both  for  .V  =  I  )  and 
0.5.  This  is  remini  cent  of  the  fact  that  in  numerous  other  studies  on  naturally  excited 
cold  air  jets.  f<)  /  for  the  most  energetic  shear  layer  mode  observed  was  likewise  some 
30-40%  below  the  predicted  value  (c.g.  Gutmark  &  ifo  (1983). 

To  understand  how  the  oscillating  mode  tends  to  dominate  the  flow  field  in  spite  of 
the  fact  that  it  is  not  associated  with  spectacular  growth  rates,  we  examine  in  figure 
24 (a)  the  strcamwisc  evolution  of  the  oscillating  mode  between  s/D  ~  1  and  b\ 
partitioning  the  interval  into  f  our  smaller  subintervals,  and  evaluating  the  amplification 
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curves  for  eat4’  of  the  subintervals.  For  each  subinterval,  the  total  amplification  of  the 
oscillating  mode  (560  Hz)  is  not  sharply  distinguished  from  the  amplification  at 
neighbouring  frequencies.  The  net  total  amplification  between  x/O  =  1  and  75  may  be 
viewed  by  comparing  the  power  spectra  obtained  at  the  two  locations,  as  seen  in  figure 
24(6).  Figures  24  (a.  h)  together  show  that  for  D/0—  83.  the  large  base-to-peak 
intensity  of  the  oscillating  mode  is  established  right  at  the  nozzle. 

Our  observations,  in  summary,  are  that,  for  D/0  —  83.  the  base-to-peak  intensity 
measured  in  the  immediate  vicinity  of  the  nozzle  is  generally  close  in  value  to  that 
measured  at  other  locations  upstream  of  the  wave-breaking  point:  quite  generally,  in 
fact,  the  oscillating  mode  is  not  a  consequence  of  an  anomalous  spatial  growth.  In  the 
next  section,  we  examine  the  effect  that  acoustic  signals  and  other  controlled  changes 
to  the  environment  have  on  the  oscillating  mode. 


6.  Response  to  acoustic  forcing  and  other  changes  in  the  environment 

6. 1 .  Sensitivity  to  coherent  disturbances 

It  has  been  established  that  the  Strouhal  number  of  the  oscillating  mode  is  independent 
of  the  maximum  initial  shear-layer  disturbance  level  u'/l  even  when  u'/i varies  by 
more  than  a  factor  of  1  (see  figure  1 3)  We  can  therefore  say  that  the  frequency  selection 
is  determined  uniquely  by  the  overall  flow  configuration,  and  by  S  and  D/0.  Figure  18 
(sj 5.2)  shows  similarly  that  the  near-field  base-to-peak  intensity  is  not  systematically 
affected  by  a'/C\.  over  the  parameter  region  of  interest.  This  means  that  a  large 
compone  it  of  the  oscillation  intensity  is  determined  by  die  overall  flow,  as  distinct 
from  external  sources.  This  appears  to  be  so  in  spite  of  the  significant  random  scatter 
seen  in  figure  18. 

The  following  experiment  is  aimed  at  helping  us  understand  the  influence  that  time- 
dependent.  spatially-coherent  disturbances  can  have  on  the  oscillation  intensity.  The 
disturbance  is  a  sinusoidal  acoustic  signal  which  is  spatially  coherent  in  the  sense  that 
the  acoustic  wavelength  is  much  larger  than  the  spatial  scales  of  the  flow :  I  D/ a ,  d  \  . 
Spatially  coherent  disturbances  could  arise  in  ‘unforced"  experiments  as  well,  either 
from  acoustic  cavity  resonances  of  the  upstream  settling  chamber,  or  from  uncontrolled 
tar-field  acoustic  sources. 

6.1.1  Overall  response 

In  figure  25.  the  r.m.s.  value  of  the  hot-wire  signal.  is  plotted  as  a  function  of  x/O 
for  different  values  of  the  acoustic  forcing  strength.  // .  p  is  actually  calculated  using 
the  voltage  which  drives  the  loudspeaker,  which  is  directly  proportional  to  the  radiated 
acoustic  pressure.  These  curves  in  figure  25  provide  only  a  qualitative  comparison  of 
fluctuation  intensities  at  differing  streamvvise  locations,  because  even  close  to  the 
nozzle,  the  mean  velocity  and  the  mean  helium  concentration  differ  for  each  x/O.  and 
so  the  proportionality  constant  relating  /:  and  a'  changes  from  place  to  place. 

Figure  25(a)  slews  the  response  of  the  heterogeneous  jet  when  forced  at  a  frequency 
(860  Hz)  that  is  unrelated  to  the  oscillating  mode  (740  Hz).  For  the  range  of  forcing 
amplitudes  considered,  the  860  Hz  mode  grows  in  a  roughly  exponent!  >1  F.shion  until 
nonlinear  damping  inhibits  further  amplification.  The  growth  rate  in  the  near  field 
( \  / 0  r:  50)  appears  to  be  unaffected  by  the  I'ororg  amplitude:  the  forcing  merely 
augments  the  spectral  intensities  in  a  spatially  uniform  manner.  The  maximum 
intensity  attained  with  dow  nstream  distance  appears  to  be  nearly  independent  of  the 
forcing  amplitude,  except  that  the  i  \imum  is  reached  closer  to  the  nozzle  for  larger 
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Fan  ki  25.  The  magnitude  of  the  signal  from  a  hot  wire  located  in  the  'hear  later,  plotted  as  a 
function  of  x  >i  for  carious  acoustic  pressure  let  els  />  m.  no  forcing,  "...  3(1  JB:  .  5(1  dB.  •  .  ~(i  d  H 
S  -  0.45;  /)■  n  -  ,S3.  (n>  Pure  tone  forcing  at  S60  II/.  and  (At  pure  tone  forcing  at  ”40  11/  (oscillating 
mode). 

forcing  amplitudes.  Acoustically  excited  shear-layer  modes  in  air  jets  exhibit  these  same 
response  characteristics  (Frymuth  1966). 

Figure  25(b)  shows  that  these  qualitative  results  also  hold  when  the  same  jet  is 
excited  at  the  frequency  of  the  oscillating  mode.  The  two  principal  differences  are :  first, 
the  intensity  of  the  oscillating  mode  in  the  absence  of  forcing  is  significantly  greater 
than  the  860  Hz  mode:  secondly,  the  fractional  increase  in  intensity  induced  b>  the 
controlled  forcing  is  correspondingly  less  for  the  oscillating  mode. 

Figure  25  helps  to  explain  previous  studies  in  acoustically  forced  heterogeneous  jets 
The  present  authors  (Kyle  16X6;  Sreenivasan  a  al.  IdSO)  found  that  at  a  fixed  location 
in  the  flow,  the  oscillating  mode  was  essentially  insensitive  to  acoustic  forcing,  while 
other  shear-layer  modes  exhibited  linear  dependence  on  /<’  when  />'  was  small  For  the 
present  conditions,  this  behaviour  is  found  at  x/0  =  100.  where  the  oscillating  mode 
response  is  nearly  independent  of  the  forcing  amplitude  (figure  25b);  whereas,  the 
shear-layer  modes,  because  they  are  initially  less  intense,  show  definite  dependence  on 
acoustic  forcing  when  the  forcing  is  small  or  moderate  (figure  25<d. 

6.1.2.  Excited  portion  of  the  response 

In  order  to  obtain  a  more  detailed,  quantitative  picture  of  the  response  to  external 
forcing,  it  is  convenient  to  decompose  the  velocity  field  as  follows: 

u'  —  »N(v.  r.  low-level  disturbances) +  nj  (  v.  r.  />  ).  (<>.l) 

Here,  it'  is  the  measured  intensity  of  velocity  fluctuations  in  the  presence  of  external 
acoustic  torcing.  tts  is  the  response  of  the  jet  with  no  externally  applied  acoustic 
forcing:  it  is  dependent  upon  the  partial  coordinates.  \  and  r.  and  possibly  upon  low- 
level  disturbances  arising  from  within  the  jet  facility  and  the  ambient,  t/,  is  the  excited 
portion  of  the  response,  and  is  dependent  upon  the  spatial  coordinates  and  upon  p  . 
This  decomposition  is  essential  for  recognizing  the  belnmour  of  u\  which  alone  is  a 
function  of  the  forcing  amplitude  whenever  us  is  of  comparable  magnitude:  this 
situation  occurs  when  the  jet  is  forced  at  the  frequency  of  the  oscillating  mode  Note 
that  (6.1)  neither  requires  nor  implies  that  receptivity  is  a  linear  process. 
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Sound  power  (dB) 

Fiot/Rf:  26.  Excited  portion  of  the  response  //,  as  a  (unction  of//:  Q.  xiO  =  1 :  .v/«  =  10; 

,\/n  —  25;  x  .  x/0  =  75;  +.  .x/0  -  200.  5  =  0.45;  D/a  =  83.  (</)  Pure  tone  forcing  at  86(1  Hz.  and  (6) 
pure  tone  forcing  at  740  Hz  (oscillating  mode).  The  arbitrary  reference  for  the  dB  scale  varies  from 
curve  to  curve. 


Figure  26(a)  shows  plotted  as  a  function  of  //  at  various  streamwise  locations  at 
860  Hz  for  the  same  flow  conditions  as  in  Figure  25.  It  is  expected  that  E'  x  a'  for  data 
obtained  at  fixed  probe  locations  very  dose  to  the  nozzle,  say  for  x/0  -  1.  where  the 
acoustically  excited  disturbances  are  small  enough  that  they  do  not  significantly  affect 
the  mean  flow.  The  data  show  that  for  excitation  levels  spanning  approximately  two 
decades,  uK  is  proportional  to  the  forcing  amplitude  at  x/0  =  1.  10  and  25.  It  is  seen 
that  the  linear  portion  of  these  three  curves  are  all  parallel  to  one  another  with  slope 
equal  to  unity.  Hence  they  must  satisfy 

log,,,  uy  =  log,,,//  +  G(.v).  (6.2) 

where  G(.v)  does  not  depend  upon  //.  Letting  G’(.v)  =  log,„g(  v).  (6.2)  may  be  rewritten 
as 

u\ •.=/>.?(  v),  (6.3) 

where  gtv)  describes  the  dependence  of  uv  on  v.  From  (6.3)  it  is  clear  that  receptivity 
is  indeed  a  linear  phenomenon  for  heterogeneous  flows,  just  as  it  is  for  homogeneous 
shear  flows.  To  our  knowledge,  this  result  has  not  been  shown  previously.  In  fact,  (6.3) 
is  a  generalization  of  the  results  obtained  for  air  jets  by  Freymuth  (1966)  who  found 
jcf(.v)  -  Be",  with  B  and  x  dependent  upon  St„  alone.  Although  the  amplification 
function  g(.v)  cannot  be  evaluated  in  this  case,  (6.3)  shows  that  growth  rate  in  the  near 
field  for  the  860  Hz  mode  is  independent  of  the  initial  disturbance  intensity.  This  result 
is  reflected  in  Figure  25(a),  where  the  amplification  curves  arc  all  parallel  to  one  another 
in  the  near  Field. 

Figure  26(b)  shows  »j.  as  a  Function  of  /?'  when  the  forcing  frequency  is  set  equal  to 
the  oscillating  mode  frequency.  Note  that  these  conditions  arc  quite  far  from  the  onset, 
as  A.V  *  0.25.  It  is  seen  that  u\,  is  approximately  proportional  to  //  at  x/0  —  I. 
furthermore,  for  these  flow  conditions,  the  data  suggest  that  (6.2)  and  (6,3)  obtain 
with  only  slight  error  over  a  region  extending  from  x/O  =  I  to  a  location  somewhere 
between  x/o  25  and  x/0  =  75.  Thus,  background  disturbances  can  affect  the 
evolution  of  vortical  structures  associated  with  the  oscillating  mode  in  much  the  same 


Instability  and  breakdown  of  a  round  variuhlc-dcnsin  jet  (At 

way  as  they  do  in  air  jets.  Indeed  il  such  disturbances  are  very  strong,  figure  25  suivjcsts 
that  they  may  also  aflect  the  wave-breaking  length  and  the  overall  structure  of  the 
mean  field  in  the  transition  region 

Finally,  we  recall  experimental  results  presented  in  Sreenivasan  et  nl.  (11>XM)  lor  the 
case  of  extremely  intense  sinusoidal  forcing.  There,  it  was  shown  that  the  Urns  can 
exhibit  behaviour  that  is  characteristic  of  a  coupled  nonlinear  oscillator,  in  that  the 
oscillatory  instability  can  become  entrained  or  Mocked-in"  to  the  frequence  of  the 
forcing. 

6.4.  Test  for  the  sensitivity  to  the  turbulence  level  at  the  nozzle  e\it 
We  have  shown  earlier  that  the  frequency  of  the  oscillating  inode  is  independent  of 
weak  ambient  disturbances.  In  the  following  experiment,  we  examine  w  hether  strong, 
spatially  incoherent  vortical  fluctuations  originating  upstream  of  the  nozzle  have  a 
strong  effect  on  the  stability  of  the  oscillating  mode  (see  also  Slrykowski  &  Russ  I W2). 

Tw  o  screens  of  differing  mesh  sizes  were  installed  separately  into  the  throat  of  the 
13.3  mm  nozzle,  and  the  resulting  qualitative  changes  in  the  power  spectrum  were 
observed  using  a  hot  wire  located  on  the  jet  centreline.  Although  the  screens  were  not 
fully  characterized,  the  turbulence  levels  along  the  centreline  in  the  nozzle  plane  were 
measured.  The  results  for  Re  =  3380  were:  (a)  with  no  screen.  //'/(,.  =  (1.003:  ih)  with 
the  small-mesh  screen.  u'/L\,  =  0.022;  (c)  with  the  large-mesh  screen,  u  .  I.  j.  ~  0.055. 

The  nozzle  fluid  was  then  changed  to  helium  while  maintaining  the  same  Re.  and  the 
three  spectra  shown  in  figure  27  were  obtained.  With  no  screen,  the  fundamental  rises 
25  dB  above  the  background  noise,  with  prominent  harmonics.  With  the  small-mcsh 
screen,  the  fundamental  rose  only  20  dB  and  is  somewhat  broadened.  With  the  large- 
mesh  screen,  the  fundamental  rose  25  dB  above  the  noise  and  retained  its  •  spikiness': 
however,  its  harmonics  were  attenuated  relative  to  the  flow  without  the  screen.  Thus, 
even  at  moderately  high  turbulence  levels  (of  the  sort  not  found  in  respectable  jet 
facilities),  the  discrete  frequency  nature  of  the  instability  was  not  quenched  For  both 
screens  the  frequency  decreased.  This  is  not  surprising  when  one  considers  that  the 
turbulent  boundary  layer  must  be  thicker  than  the  laminar  one.  and  so  by  the  scaling 
laws  described  in  §5.  the  frequency  is  expected  to  decrease. 

6.5.  Test  for  the  sensitivity  to  the  ambient  environment 
Another  experiment  was  made  to  test  the  sensitivity  of  the  oscillating  mode  to  the 
spatial  confinement  of  the  jet.  If  the  jet  is  surrounded  by  an  open-ended  circular 
container  sitting  on  the  nozzle  block,  the  oscillating  mode  hardly  changes  in  intensity 
but  shows  a  frequency  increase  by  about  6  or  7%  when  the  container  diameter  and 
height  are  about  ten  nozzle  diameters.  When  the  top  of  the  container  is  covered  In  a 
flat  board  with  a  small  central  hole  (diameter  on  the  order  of  the  nozzle  diameter),  the 
oscillating  mode  disappears;  this  is  not  surprising  because  the  recirculating  flow  cot  up 
in  the  container  renders,  in  due  course,  the  density  of  the  ambient  gas  equal  to  that  of 
the  nozzle  gas.  Experiments  with  several  geometrical  combinations  of  the  container 
have  suggested  to  us  that  the  oscillating  mode  disappears  only  under  drastic  changes 
of  the  environment. 

6.6.  Sensitivity  to  local  perturbations 

While  the  oscillating  mode  is  robust  to  many  types  of  perturbations,  it  is  quite  sensitive 
to  modifications  of  the  mean  velocity  field.  A  case  in  point  is  its  behaviour  in  the 
presence  of  an  external  body  such  as  a  small-diameter  pin  located  at  a  suitable  place 
inside  the  flow.  This  was  investigated  by  Sreenivasan  et  at.  (!5>8‘>)  who  placed  in  the 
centre  of  the  jet  perpendicular  to  the  axis  a  straight  pin  of  a  certain  diameter  and 


654  D.  M .  Kyle  and  K.  R.  Sreenintsun 


o  2000 

i 

Frequency  (H/) 

Fkh  ri.  27.  Power  spectra  showing  the  influence  of  grid-generated  turbulence  on  the  oscillating 
mode:  .S'  =  0.14:  D/ft  =  55:  x/D  =  2:  r  =  0:  D  —  13.3  mm.  (</)  No  grid,  (/'/(  =  0.003 ;  (b)  finer  grid. 
tt'/V,  -  0.022:  (o)  coarser  grid,  if'/U,.  =  0.055. 

monitored  the  oscillating  mode  as  a  function  of  the  pin  position  along  the  axts. 
Typically,  the  pin  diameter  was  an  order  of  magnitude  smaller  than  the  nozzle 
diameter.  They  observed  that  the  oscillating  mode  could  be  suppressed  whenever  the 
pin  was  located  in  a  certain  neighbourhood  in  the  flow.  The  observation  suggests  the 
possibility  that  the  oscillating  mode  is  a  consequence  of  a  local  instability  of  the  flow 
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in  that  neighbourhood,  and  that  it  disappear  because  the  pin  significantly  disrupts  the 
mean  flow  near  where  the  instability  arises  (see  also  Monkevvitz  et  al.  1990).  It  should 
be  noted  that  perturbations  to  the  mean  velocity  profile  can  affect  the  jet  development 
in  isothermal  air  jets  as  well  (Bradbury  &  KhuJem  1975). 


7.  Summary  and  discussion 

7.1.  Distinction  between  various  modes 

In  heterogeneous  jets  which  do  not  support  the  oscillating  mode,  shear-layer 
disturbances  evolve  in  an  analogous  manner  to  those  in  constant  density  jets.  We  have 
found  that  very  close  to  the  nozzle,  the  power  spectrum  of  u  is  determined  by  a 
combination  of  the  background  spectral  peaks  at  the  nozzle  exit  and  the  spatial 
amplification  rate  of  shear-layer  modes.  The  latter  is  a  smoothly  varying  function  of 
the  frequency  of  the  disturbance.  The  maximum  of  this  function  is  the  'most  spatially 
amplified  mode’.  The  most  spatially  amplified  mode  is  a  smoothly  varying  function  of 
S.  The  non-dimensional  frequency  JO/  i/,  of  this  mode  is  in  agreement  with  other 
experimental  studies  in  air  jets,  fd/ly  decreases  with  .S',  as  predicted  by  the  spatial 
stability  theory  for  non-uniform  density  jets.  From  films  we  found  that  the  average 
value  of  the  wavelength  A/0  is  also  well-predicted  by  the  theory. 

Important  changes  occur  whenever  the  jet  supports  the  oscillating  mode.  The 
structure  in  the  near-field  is  highly  organized  both  spatially  and  temporally,  and 
centreline  values  of  u'/b\,  can  become  extraordinarily  large.  The  processes  of  wave 
breaking  and  roll-up  formation  occur  relatively  close  to  the  nozzle.  The  non- 
dimensional  intensity  of  the  oscillating  mode  and  Strouhal  number  are  definite 
functions  of  S  and  D/0  alone,  that  is  they  are  not  systematically  affected  by 
background  disturbances  or  by  Re  within  the  parameter  ranges  of  this  study.  The 
‘most  spatially  amplified'  shear-layer  mode  also  depends  smoothly  on  S  and  D/0,  but 
the  characteristics  {A/0.f0j (./,)  of  this  mode  have  been  shown  to  be  unrelated  to  those 
of  the  oscillating  mode.  The  base-to-peak  intensity  of  the  oscillating  mode  measured 
in  the  immediate  vicinity  of  the  nozzle  is  generally  close  in  value  to  measurements  made 
further  downstream  in  the  near  field.  Thus,  the  shear  layer  acts  as  a  spatial  amplifier 
in  the  usual  way;  the  oscillating  mode  does  not  exhibit  an  anomalously  large  spatial 
growth  rate.  The  initial  intensity  of  the  oscillating  mode  can  be  influenced  by  acoustic- 
forcing  in  a  linear  fashion,  but  in  a  quiet  laboratory  this  component  is  small  relative 
to  the  flow-induced  portion. 

7.2.  Vort icily 

Although  we  have  not  measured  the  vorticity  distribution  in  the  presence  of  the 
oscillating  mode,  it  is  possible  to  argue  that  the  oscillating  mode  is  associated  with 
highly  concentrated  vortical  structures.  The  instability  of  small  disturbances  in  a 
laminar  jet  results  in  the  migration  of  vorticity  to  form  periodic  concentrations.  In  a 
qualitative  sense,  the  flow  field  which  accompanies  such  vorticity  disturbances  is 
analogous  to  the  nonlinear  solutions  of  the  vorticity  equation  introduced  by  Stuart 
(1967).  These  solutions,  which  are  periodic  in  one  direction  and  have  shear  in  the  other, 
were  used  by  Stuart  to  describe  the  flow  patterns  in  a  shear  layer  with  periodic  vorticitv. 
The  stream  function  »//  for  these  solutions  is  given  by 

i//  ~  er  +  In  [Ccosh  (r)  +  A  cos  (  v  —  <■/)].  (7.1) 

where  c  is  the  wave  speed,  and  A  and  (’  are  related  bv 

A  ~  {C2-  1)’\ 


(7.2) 


656 


D.  M.  Kyle  and  K.  R.  Sreenirasan 


The  parameter  C  indicates  the  degree  of  concentration  of  the  vorticity:  ('=  1 
corresponds  to  uniform  vorticity  in  the  \-direction.  while  (  ■  x  corresponds  to  a  row 
of  point  vortices  (  Lamb  1945).  Stuart  calculated  u'/U,.  using  (7.1)  and  (7.2)  and  found 
that  at  every  point  in  the  flow,  the  intensity  increases  monotonically  with  Note  that 
the  area  integral  of  vorticity  over  any  fixed  region  does  not  change  with  C  in  Stuart's 
solution:  only  the  distribution  is  affected.  The  conclusion  is  that  for  a  given  amount  of 
mean  vorticity.  the  velocity  fluctuation  intensity  can  be  strongly  affected  by  the  degree 
to  which  the  vorticity  disturbance  is  concentrated.  We  have  measured  the  r.m.s. 
velocity  u‘ / L\,  along  the  centreline  for  two  jets  with  density  ratios  S  =  1 .0  and  S  -  0.29 
but  with  identical  (non-dimensional)  mean  vorticity  distribution  at  the  nozzle  exit 
(vj4.4).  We  found  that  the  jet  supporting  the  oscillating  mode  (.V  =  0.29)  exhibits 
significantly  larger  centreline  fluctuation  intensities  (figure  12).  This  suggests  that  the 
vorticity  is  much  more  concentrated  in  the  streamwise  direction  in  jets  w  hich  support 
the  oscillating  mode. 

Winant  &  Brow'and  (1974)  have  proposed  a  model  for  the  vortex  pairing  process 
based  on  Stuart's  solution.  Their  data  show  that  for  a  given  wave  speed,  both  the  initial 
growth  rate  of  the  subharmonic  instability  and  the  co-rotating  speed  of  pairing  vortices 
will  increase  with  C.  If  it  is  true  that  the  vorticity  is  more  highly  concentrated  for  the 
oscillating  mode,  then  Winant  &  Browand's  results  would  suggest  that  for  S  =  0.29. 
the  pairing  process  should  occur  sooner  and  that  the  process  should  be  more  energetic. 
We  have  indeed  shown  (see  table  1)  that  the  pairing  generally  occurs  sooner  in  jets 
dominated  by  the  oscillating  mode,  and  that  the  fluctuation  levels  reaches  abnormally 
intense  values  (figure  12)  near  the  location  of  vortex  pairing. 

7.3.  Theoretical  models 

7.3.1.  Linear  theory 

We  have  not  succeeded  in  experimentally  identifying  the  direct  physical  cause  of  the 
unusual  instability  arising  under  certain  circumstances  in  variable-density  jets.  Various 
measurements  discussed  in  §§4  and  5  show  that  a  plausible  explanation  is  provided  by 
the  linear  theory  for  wave  packet  growth.  Agreement  and  disparity  between  theory  and 
experiment  will  now'  be  reviewed. 

If  oy(k*(0))  >  0,  the  theory  shows  that  in  a  streamwise  homogeneous  system,  any 
spatially  localized  disturbance  will,  as  time  increases,  lead  to  an  increasingly 
monochromatic  response,  as  observed.  In  fact,  for  S  —  0.50  there  is  a  reasonable 
correspondence  between  the  predicted  and  observed  values  o ffO/L/v,  within  a  limited 
range  of  D/0  values  around  80.  In  this  same  range  of  D/O,  approximate  correspondence 
has  been  found  between  the  predicted  critical  value  of  ,S’(=  0.72)  and  the  observed 
onset  value  S„  ( =  0.61). 

On  the  other  hand,  as  D/0  increases  much  above  about  80.  the  correspondence  with 
the  theory  breaks  down.  In  contradiction  to  the  theory,  S„  rapidly  takes  smaller  values 
with  increasing  D/0  (figure  19).  For  large  enough  values  of  D/0  the  oscillating  mode 
vanishes,  even  though  the  theory  predicts  that  for  all  ,S  -<  0.62.  w,(  A  *(()))  >  0  even  in  the 
limit  of  0  -0  (Monkewitz  &  Sohn  1986,  1988).  A  complete  validation  of  any  linear 
stability  theory  might  comprise  a  matching  of  velocity  data  with  computed 
eigenfunctions  and  eigenvalues,  and  a  matching  of  any  predicted  critical  parameter 
values  with  their  experimentally  determined  onset  values.  By  these  standards,  the 
correspondence  between  the  oscillating  mode  and  the  spatio-temporal  theory 
summarized  in  §  I  would  seem  quite  scant. 

Even  the  somewhat  limited  correspondence  found  in  and  5  between  the  linear 
theory  and  experiment  may  seem  surprising  because:  (</)  the  instability  theory  assumes 
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homogeneity  in  the  streamwise  direction,  whereas  the  How  in  the  near  field  is  strongly 
inhomogeneous,  and  (b)  the  oscillating  mode  is  so  intense  that  linear  mechanisms  can 
be  called  into  question  without  much  thought.  As  regards  (a),  we  have  shown  that  the 
oscillating  mode  is  global  in  the  sense  that  the  separate  processes  involved  in 
laminar  turbulent  transition  occur  in  a  globally  organized  fashion.  On  the  other  hand, 
its  initiation  may  be  related  to  the  local  spatio  temporal  instability  of  a  certain  profile 
in  the  near  field.  Some  support  for  this  for  this  conjecture  comes  from  Chomaz.  Huerre 
&  Redekopp  (1988)  who  have  shown  that  it  is  possible  to  construct  a  one-dimensional 
linear  system  which  exhibits  just  this  behaviour.  They  have  shown  that  if  a  certain 
coefficient  of  the  linearized  Ginzburg  Landau  equation  is  allowed  to  vary  spatially, 
then  within  the  context  of  that  system,  the  existenee  of  temporally  growing  modes  of 
the  form  <;5(.v) ei'"t  always  requires  that  m,(A*(0))  >  0  for  some  spatial  interval  (see  the 
review  by  Huerre  &  Monkewitz  1990).  Concerning  (h),  a  common  experience  is  that  the 
frequency  of  instability  between  the  linear  and  nonlinear  stages  does  not  vary  much. 
We  stress  that,  although  the  linear  theory  might  help  in  understanding  the  behaviour 
of  the  instability,  large-amplitude  motions  play  an  essential  role  in  sustaining  the  stable 
periodic  state  overall. 

7.3.2.  Nonlinear  theory 

Raghu  &  Monkewitz  (1991)  have  pointed  out  that  considerable  physical  insight 
might  be  obtained  from  studying  the  temporal  growth  of  the  oscillating  mode.  They 
have  tried  to  show  that  the  oscillatory  disturbance  field  in  a  heated  air  jet  can  be 
modelled  as  a  ‘global  temporal  mode',  that  is. 

u'(.\\t)  =  A(l)f(x).  (7.3) 

and  that  this  mode  obeys  the  Landau  equation.  Because  of  the  similarity  between  the 
two  flows,  such  a  conclusion  would  also  be  appropriate  to  the  heterogeneous  flows 
considered  here.  Unfortunately,  their  results  are  inconclusive  in  that  they  obtained 
data  at  only  a  single  streamwise  location  and  provide  no  information  as  to  whether  the 
factorization  (7.3)  is  valid.  Moreover,  our  examination  of  onset  for  a  wide  range  of 
flow  conditions  (§5.4)  indicates  that  the  spatial  distribution  of  fluctuation  intensity  can 
itself  be  strongly  amplitude-dependent  near  onset.  For  such  behaviour  the  factorization 
(7.3)  is  inappropriate.  Yet,  the  ‘bifurcation'  appears  to  be  of  the  supercritical  type  for 
jets  produced  in  many  laboratory  facilities.  (For  an  exception  under  certain 
circumstances,  see  Sreenivasan  et  al.  1989.)  In  further  pursuit  of  this  line  of  inquiry,  we 
propose  that  transient  response  simultaneously  measured  at  different  streamwise 
locations  would  be  very  valuable.  At  least  under  flow  conditions  with  slow  transient 
growth,  one  might  learn  more  about  the  incipient  stages  of  the  oscillating  mode. 

We  thank  Professors  Edward  Bolton  and  Boa-Tch  Chu  for  useful  discussions.  The 
work  was  financially  supported  by  a  grant  from  the  Air  Force  of  Scientific  Research. 


Appendix.  Interpreting  hot-wire  spectra  of  small  disturbance  in  laminar 
heterogeneous  shear  layers 

While  one  is  primarily  interested  in  measuring  the  power  spectral  density  of  a'  in  the 
shear-layer  region,  a  single  hot  wire  cannot  be  used  in  regions  of  variable  density  w  here 
the  mass  fraction  of  helium  c  fluctuates.  This  is  because  the  hot-wire  voltage  I.  is 
affected  by  both  fields: 


/:( I :  v.  /•)  =  /•.'( r( t ;  v.  /•).  U(t:  v.  /  ) ). 


(A  I) 
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Frequency  (Hz) 

Figure  28.  Power  spectra  of  u'  obtained  at  different  radial  locations  in  the  shear  layer  region:  S  = 

1.00:  DjH  =  83:  x/0  -  15: - .  potential  core: -  .  r  =  \D  +  50.  The  arbitrary  reference  for  the 

dB  scale  varies  for  the  two  curves. 


Instead,  the  power  spectral  density  of  u'  has  been  obtained  using  a  hot  wire  placed 
within  the  ‘potential  core'  region,  where  both  the  mean  velocity  and  the  mean  helium 
concentration  are  uniformly  distributed  in  the  radial  direction.  Because  c  can  only 
decrease  owing  to  mixing,  it  follows  that  concentration  fluctuations  are  altogether 
absent  in  the  potential  core.  For  each  of  the  flows  discussed  in  this  Appendix,  the 
existence  of  this  region  has  been  confirmed.  In  each  example,  the  hot  wire  was  placed 
very  near  to  where  either  U  or  c  (or  both)  is  just  beginning  to  decrease  with  r.  i.e.  near 
the  ‘inner  edge’  of  the  shear  layer.  At  this  location  the  hot  wire  responds  only  to  the 
velocity  fluctuations  associated  with  shear-layer  disturbances,  which  are  presumed  to 
be  strong  along  the  inner  edge  of  the  shear  layer  (see  Freymuth  1966). 

Flow  uniformity  within  the  potential  core,  as  well  as  the  location  of  the  inner  edge 
of  the  shear  layer,  was  confirmed  in  each  case  as  follows.  As  the  hot  wire  is  moved 
radially  outward  from  the  centreline  by  small  increments  Hr.  the  mean  voltage  changes 
according  to: 


E(r  +  Ar)  —  E{r)  = 


<EcU  cEtc 

(f^  + -  (,)>). 
<  {/(  /•  it  tr 


(A  2) 


The  partial  derivatives  <  E/CU  and  eE/Cc  have  been  measured  directly  by  placing  the 
hot  wire  at  the  centre  of  the  nozzle  exit  and  separately  varying  U  and  c  for  several 
different  conditions.  It  was  found  that  I'E/eU.  eE/ec  >  0.  while  (CE/H U)/(eE/ee)  is  of 
order  unity  (see  also  Way  &  Libby  1971 ).  Because  ('(.//<>,  <Y/fr  <  0  the  mean  voltage 
will  always  decrease  with  r.  Thus,  the  radial  extent  of  the  potential  core  can  be 
unambiguously  determined. 

Unfortunately,  it  is  not  convenient  to  rely  upon  measurements  obtained  in  the 
potential  core  because  the  presence  of  the  probe  itself  can  alter  the  evolution  of  shear- 
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f'KiURi;  29.  Power  spectra  of  u  obtained  at  different  radial  locations  in  the  shear  layer  region:  .V  = 

0.70:  D/0  ~  83: - .  potential  core:  ,  r  =  \D  +  5f)  («)  x/0  =  25;  (b)  x/t)  =  75.  The  arbitrary 

reference  for  the  dB  scale  varies  for  the  two  curves. 


Instability  and  breakdown  of  a  round  variable-density  jet 


66 1 


Fiona.  31.  Mode  amplification  between  x/tt  --  25  and  x/tt  —  75  as  a  function  of  frequency  for  two 
different  radial  locations:  S  —  0.60.  D/It  =  62.5:  .  potential  core:  .  r  =  ID  +■  5"  The  gain 

at  each  r  is  determined  to  within  an  arbitrary  constant  factor. 

layer  disturbances  (Hussain  &  Zatnan  1078;  Sreenivasan  el  a/.  1989).  We  have 
therefore  measured  the  spectrum  near  the  outer  edge  of  the  shear  layer  at  r  =  !,D  +  5d 
and  exploited  the  result,  discussed  below,  that  in  regions  near  the  nozzle  where 
disturbances  are  small,  the  power  spectral  density  function  obtained  at  r  =  \D  +  5b  is 
approximately  the  same  as  the  power  spectrum  of  if  obtained  in  the  potential  core, 
except  for  a  constant  factor  which  is  independent  of  the  frequency. 

Figure  28  shows  power  spectra  of  it'  obtained  in  an  air  jet  within  the  potential  core 
(r  —  -,D—  10(7)  and  near  the  outer  edge  of  the  shear  layer  (r  =  \D  +  5d)  at  x/0  =  75.  The 
spectra  have  been  multiplied  by  arbitrary  constants  selected  to  show  the  similarity  in 
the  form  of  th.  curve.  It  is  seen  that  the  spectra  are  very  similar,  except  at  the  lowest 
frequencies.  Measurements  in  air  jets  obtained  at  different  streamwise  locations  x/0  < 
75  and  for  differing  D/d  all  yield  similar  results.  Figure  29  shows  power  spectra  for 
S  =  0.70  at  v  =  25(7  and  75(7.  Again  the  spectra  obtained  within  the  potential  core  and 
at  r  =  if) +  5(7  are  approximately  the  same  except  for  a  constant  normalizing  factor. 
Figure  30(a)  shows  that  this  similarity  can  extend  even  to  jets  which  support  the 
oscillating  mode,  so  long  as  measurements  are  made  close  to  the  nozzle.  In  figure  30(6). 
the  similarity  has  degraded  somewhat,  presumably  owing  to  large-amplitude  effects. 
Note  that  for  each  potential  core  measurement,  possible  effects  associated  with  the 
intruding  probe  were  monitored  using  a  second  hot  wire  located  well  outside  the  shear 
layer.  No  intrusion  effects  were  observed  for  the  examples  shewn  in  figures  28  30. 

During  the  course  of  this  study,  spectra  obtained  at  r  =  \D  +  5d  are  used  for 
calculating  the  following  measures:  (i)  base-to-peak  intensity  measured  in  the 
immediate  vicinity  of  the  nozzle  (x/D  —  0.10).  Whenever  probe  intrusion  is  not  a 
factor,  this  measure  is  found  to  be  roughly  independent  of  r  (see  figure  30c/).  (ii)  Most 
spatially  amplified  mode  (§4.4).  Figure  29  suggests  that  the  functional  relation  between 
spatial  amplification  rate  and  frequency  should  be  approximately  independent  of  the 
radial  position,  except  for  a  constant  factor  that  is  independent  of  frequency.  More 
particularly,  the  maximum  of  the  spatial  amplification  curve  (§5.5)  should  be 
independent  of  radial  position.  This  is  confirmed  in  figure  31  for  S  =  0.60  and  D/d  = 
62.5.  Again,  the  spectra  arc  seen  to  differ  most  in  form  at  the  lowest  frequencies,  (iii) 
Frequency  of  the  oscillating  mode.  This  measure  is  actually  a  constant  for  all  points  in 
the  transitional  flow. 
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Abstract — -The  OH  radicals  in  a  hydrogen  diffusion  flame  burning  in  ambient  air  wm*  sisyah/ed  using 
planar  laser-induced  fluorescence.  Local  thicknesses  of  the  OH  regions  were  measured  at  three  downstream 
distances.  Measurements  show  that  the  thickness  is  distributed  approximately  according  to  a  lognormal 
probability  density  function. 


1.  INTRODUCTION 

The  spatial  structure  of  turbulence  plays  u  significant  role  in  determining  dynamical 
processes  such  as  mixing,  stability,  and  lift-off  in  turbulent  flames.  Understanding  the 
basic  structure  of  turbulent  diffusion  flames  involves  many  issues  such  as  turbulent  mixing, 
scalar  diffusion,  multistep  chemical  kinetics  and  heat  release.  Molecular  mixing  occurs 
across  scalar  interfaces  between  the  fuel  and  the  oxidant.  Combustion  occurs  when  mixing 
is  completed  to  stoichiometric  limits  within  these  interfaces.  The  corresponding  heat 
release  will  in  turn  have  a  large  effect  on  the  local  structure  of  turbulence.  This  interplay 
is  quite  complex,  and  we  shall  not  be  concerned  here  with  those  details.  Instead,  we 
concentrate  on  one  aspect,  namely  the  thickness  distribution  of  hydroxyl  radical  (OH)  in 
a  hydrogen-air  diffusion  flame.  OH  is  a  product  of  combustion  and.  with  some  limitations 
to  be  discussed  later,  is  considered  to  be  a  diagnostic  of  the  flame  zone. 

The  spatial  structure  of  OH  regions  has  previously  been  characterized  in  terms  of 
thickness,  axial  scale,  orientation,  and  location  in  hydrogen-air  diffusion  flames  at  various 
Reynolds  numbers  (Seitzman  et  al.,  1990).  However,  the  probability  density  function 
representing  the  thickness  of  OH  regions  has  not  previously  been  studied  in  detail.  In 
cold  flows,  the  thickness  of  mixing  interfaces  is  on  the  order  of  the  smallest  scales  of 
motion  (LaRue  &  Libby,  1974;  Sreenivasan  et  al„  1989).  It  can  further  be  shown  that 
the  probability  distribution  of  such  scales  is  lognormal.  If  reaction  zone  structures  have 
similar  properties,  it  is  possible  that  our  knowledge  of  the  cold  flow  structure  may  be 
translated,  even  if  in  a  somewhat  limited  sense,  to  the  characterization  of  combusting 
systems. 

With  this  in  mind,  we  have  made  measurements  of  the  thickness  distribution  of  the 
OH  regions  in  a  low-to-moderatc  Reynolds  number  hydrogen  diffusion  flame  burning 
in  air.  The  measurements  show  that  the  OH  regions  are  quite  fat  (in  comparison  with 
the  Kolmogorov  thickness  of  the  cold  flow),  and  dint  their  thickness  is  approximately 
lognormal  in  distribution.  Some  remarks  are  made  on  the  interpretation  of  these 
observations. 

2.  EXPERIMENTAL  PROCEDURE 

Experiments  to  visualize  a  diffusion  flame  were  designed  at  the  United  Technologies 
Research  center.  The  goal  was  to  identify  the  reaction  region  and  obtain  a  data  base  for  its 
thickness  distribution.  Because  of  the  difficulties  in  imaging  instantaneous  reaction  zones 
of  hydrogen  flames  and  the  good  signal  qualities  of  OH.  images  of  OH  regions — which 
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had  a  capacity  of  512  <  512  pixels  and  the  overall  rest  *Ui  turn  was  .ih«ut  <J  i  mm  pixel.  Laser  sheet  thickness 
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mark  areas  where  combustion  has  been  completed — were  recorded.  The  a;  pucahility  of 
OH  as  a  marker  of  reaction  regions  in  flames  has  been  discussed  by  Hanson  ( 19N7)  and 
Vandsburger  et  al.  (1988),  and  is  also  addressed  later  in  this  paper. 

The  experimental  set  up  is  shown  in  Figure  I.  Hydrogen  was  ejected  upwards  from 
a  2.5  mm  diameter  round  nozzle  into  ambient  air  with  a  nozzle  velocity  of  108  m  s. 
The  cold  flow  Reynolds  number  based  on  nozzle  diameter  was  Rc  ~  2400.  Images  were 
recorded  at  axial  positions  centered  around  distances  of  38.  58  and  Os  nozzle  diameters. 
Some  details  are  given  below. 

The  OH  free  radicals  from  the  products  of  hydrogen-air  reaction  were  visualized  bv 
planar  laser-induced  fluorescence  (PLIF).  A  Speetra-Physics  Nd.YAG  laser  frequency- 
doubled  to  a  wavelength  of  532  nm  pumped  a  dye  laser  which  w  as  frequency -doubled  to 
a  peak  absorption  wavelength  of  OH  atoms  at  282  nm.  The  dye  laser  beam  was  passed 
through  an  optical  arrangement  which  produced  a  laser  sheet  of  thickness  of  the  order 
of  100  pm.  The  laser  pulse  width  was  on  the  order  of  10  nsec  and  the  laser  energy  at 
the  flame  was  on  the  order  of  20-30  m.l.  Fluorescence  of  the  OH  radicals  were  induced 
by  the  energy  and  wavelength  of  the  laser  sheet  where  the  (1.0)  band  of  the  A  A  '  -YU. 
transition  was  excited  at  282  nm.  This  gives  fluorescence  from  both  (1.1)  and  (0.0)  bands. 
Temperature  and  quenching  effects  on  OH  fluorcscerce  were  minimized  by  adjusting  the 
laser  to  saturate  the  transition.  Pfeflerle  et  al.  (1988)  have  demonstrated  the  effectiveness 
of  such  a  technique  in  the  PLIF  images  of  OH  in  a  chemically  reacting  boundary  layer. 
In  their  work,  saturation  was  achieved  to  within  10G.  Illuminated  regions  were  recorded 
on  an  intensified  XYBION  CID  camera  with  a  gating  synchronized  with  the  pulsing  of 
the  laser.  For  those  measurements,  the  gating  time  was  less  than  I  ps  and  the  camera 
was  operated  in  a  non-interlaced  mode.  Data  were  recorded  through  an  interference 
filter  with  a  full  width  half-maximum  (FWHM)  of  10  nm  centered  at  308  nm.  allowing 
detection  primarily  from  the  (0,0)  band.  The  camera  resolution  of  512  x  512  pixels 
allowed  an  average  image  resolution  of  0.09  mm  pixel.  Output  of  the  camera  was  fed  to 
a  video  recorder  for  continuous  acquisition. 

The  videotaped  data  were  digitized  on  a  Macintosh  11  using  a  model  D 12255  Data 
Translation  digitizing  card  with  a  640  x  480  square  pixel  format  and  associated  image 
processing  software,  images  were  digitized  to  8-bits  and  the  dynamic  range  was  on  the 
order  of  30;  however,  the  precise  value  of  the  signal  to  noise  ratio  varied  from  pixel  to 
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pixel,  being  the  lowest  where  the  structures  were  diffuse  and  OH  concentrations  low. 
The  measurements  here  were  limited  to  assessing  the  relative  thickness  of  the  spatial 
distribution  of  the  OH  regions,  and  the  absolute  value  of  pixel  intensities  and  their 
relation  to  OH  concentration  were  used  only  to  determine  the  'edges'  of  OH  regions 
Considering  the  system  as  a  whole,  the  projected  resolution  amounts  to  approximately 
0.095  mm  pixel  in  the  vertical  direction  and  0.15  mm  pixel  in  the  horizontal  direction 
over  a  46  mm  x  46  mm  area.  The  smallest  structure  measured  was  on  the  order  of  !  mm 
in  width. 

Twenty  randomly  selected  images  of  the  flame  at  each  of  the  three  streamw  ise  locations 
were  digitized.  The  thickness  of  the  OH  region  in  each  of  these  images  was  measured  at 
equal  intervals  along  the  length  of  the  flame.  At  each  interval,  the  shortest  line  across  the 
OH  region  was  defined  and  a  profile  of  the  pixel  intensities  along  the  line  was  obtained. 
The  edge  of  the  OH  region  was  identified  by  setting  a  threshold  of  pixel  intensity  from 
which  best  defined  the  width  of  the  profile  and  was  sufficiently  above  the  noise  level  of 
the  image.  The  number  of  pixels  along  the  line  crossing  the  OH  regions  with  intensities 
above  the  threshold  value  were  then  counted.  This  procedure  was  repeated  using  three 
different  pixel  intensity  thresholds,  corresponding  to  intensities  60.  70  and  SO  percent 
below  the  maximum  intensity.  From  visual  observations,  it  was  found  that  much  lower 
thresholds  would  have  been  contaminated  by  noise.  For  comparison,  hand  measurements 
were  also  taken  directly  from  images  projected  on  to  a  large  screen.  In  this  technique 
the  edges  of  the  OH  regions  were  determined  visually  rather  than  by  pixel  thresholding. 
It  should  be  noted  that  glancing  cuts  of  the  laser  through  the  jet  add  uncertainties  to 
the  interpretation  of  the  data.  Structures  from  the  images  which  were  clearly  deemed 
as  glancing  cuts  were  not  used  in  the  analysis.  From  the  recorded  data,  histograms  of 
the  thickness  of  the  OH  regions  were  produced. 

3.  RESULTS  AND  DISCUSSION 

I  he  basic  parameters  of  the  cold  flow,  estimated  according  to  the  'standard'  evolution 
of  round  jets  (e.g.,  Wygnanski  &  Fiedler.  1968).  are  shown  in  Table  1.  These  data  are 
not  directly  relevant  to  the  local  flame  front  because  the  latter  is  nearly  laniinarized  bv 
the  large  heat  release.  Yet,  the  local  reaction  .egions  are  buffeted  by  the  turbulence 
in  the  core  of  the  flow;  some  idea  of  the  cold  flow  data  was  therefore  thought  to  be 
appropriate.  The  turbulent  Reynolds  number  Re,  is  based  on  the  estimated  root-mean- 
square  (rms)  centerline  velocity  and  integral  scales.  The  Kolmogorov  length  and  velocity 
scales,  estimated  according  to  the  standard  turbulence  scaling  relationship,  arc  included, 
so  are  estimates  of  Damkbhler  number.  Da.  based  on  chemical  reaction  time  scales  of 
both  the  fon  ard  and  recombination  reactions. 
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FIGURE  2.  Visualizations  of  OH  concentration  in  a  hydrogen  diffusion  flame  at  three  locations  downstream 
of  the  nozzle.  Reynolds  number  at  the  nozzle  is  Re  =  2400.  Nozzle  diameter,  d  *  2.5  mm.  Width  of  the 
images  is  x  ~  46  mm  in  real  space. 


Before  we  present  the  results,  a  comment  is  in  order  on  the  effect  of  buoyancy,  which 
might  be  expected  to  have  increasing  influence  with  increasing  distances  from  the  nozzle. 
Some  symptoms  of  buoyancy  in  jet  flows  include  increases  in  momentum  flux.  Reynolds 
number,  and  activity  and  size  of  large  eddies.  In  studies  of  buoyancy  effects  on  turbulent 
diffusion  flames  of  propane,  Becker  &  Yamazaki  (1978)  found  that  transition  from  a 
momentum  driven  to  a  buoyant  driven  state  begins  when  a  nondimensional  streamwise 
coordinate,  £  =  Ri]/3(x/Ds),  approaches  a  value  between  1  and  2.  Here,  Ri,  is  a  'source' 
Richardson  number  based  on  an  effective  'source  diameter’  given  by  D, .  The  value  of 
£  is  substantially  smaller  in  our  flame  except  at  the  last  measuring  station,  where  it  is  of 
the  order  of  3.  Thus,  a  priori  considerations  suggest  that  some  effect  of  buoyancy  may¬ 
be  present  at  the  last  measuring  station.  We  shall  later  take  a  closer  look  at  this  issue. 

Figure  2  shows  three  typical  realizations  of  the  OH  regions  at  the  three  downstream 
positions  mentioned  earlier.  As  expected,  they  become  more  convoluted  with  increasing 
downstream  distance.  The  OH  regions  are  thin  in  comparison  with  the  flow  width,  hut 
their  thickness  is  far  larger  than  the  Kolmogorov  scale  of  the  cold  flow.  As  observed 
in  a  similar  experiment  on  a  hydrogen-air  diffusion  flame  (Kychakoff  ct  al.,  1984).  this 
relative  fatness  of  the  OH  regions  is  in  part  due  to  the  low  Reynolds  number  of  the 
flow.  In  the  experiment  by  Kychakoff  et  al.  (1984)  it  was  found  that  with  increasing 
Reynolds  number,  the  width  of  the  structures  became  smaller  and  had  greater  variations 
in  size.  Similar  characteristics  were  also  observed  by  Seitzman  et  al.  (1990).  who  note 
that  the  diffuse  structures  are  predicted  well  by  the  quasi-equilibrium  distributed  reaction 
(QEDR)  model  developed  by  Bilger  (1988).  A  detailed  review  of  the  influence  of  other 
factors  such  as  dissipation,  mixture  fraction  and  their  role  in  the  structure  and  thickness 
of  turbulent  nonpremixed  flames  can  be  found  in  Bilger  (1988). 

Because  of  the  relatively  slow  recombination  time  of  OH.  we  also  expect  the  thickness 
of  the  OH  regions  to  be  broader  than  the  instantaneous  reaction  zones.  Bilger  (1989a) 
gives  an  estimate  for  the  flame  thickness  A  in  terms  of  the  Kolmogorov  scale  >/  for  one- 
step,  second  order  irreversible  reactions  as  A/?;  =  Da~,^Rcl'4Sc~t'2'yt~h''  where  Sc  is 
the  Schmidt  number  and  ■yi  is  the  rr.is  of  mixture  fraction.  This  relationship  qualitatively 
explains  how  the  Damkohlcr  number  can  affect  the  size  of  a  reaction  zone.  As  shown  in 
Table  I,  for  hydrogen-air  reactions  Da  is  much  smaller  for  the  reverse  reaction  involving 
the  free  radical  OH  than  for  the  forward  reaction.  Hence  we  expect  OH  regions  to  cover 
a  larger  area  than  the  instantaneous  flame  region.  For  laminar  counterflow  methane-air 
diffusion  flames,  predicted  number  densities  of  both  CH  and  OH  in  mixture  fraction 
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FIGURE  3.  Distributions  of  the  thickness  of  the  OH  regions  at  38.  58,  and  98  diameters  downstream  of 
the  nozzle,  (a)-(c):  Edges  of  OH  regions  determined  by  setting  a  threshold  on  the  pixel  intensity  (see  text 
for  details),  (d)-(f):  Edges  of  OH  regions  determined  visually.  Smooth  curves  denote  fits  using  mean  and 
rms  values  from  the  measurements. 


space  (Bilger,  1989b)  have  shown  that  the  OH  regions  are  broader  and  exist  in  high 
concentrations  slightly  to  the  lean  side  of  the  instantaneous  reaction  zone  marked  by 
CH. 

Figures  3(a)-3(f)  are  histograms  of  the  thickness  of  OH  regions  taken  from  the  digitized 
images  at  38, 58  and  98  diameters  downstream.  Figures  3(a)-3(c)  represent  measurements 
obtained  by  defining  the  ‘edges’  of  the  OH  regions  by  a  threshold  which  is  60%  below 
the  largest  intensity.  It  was  found  that  the  histograms  were  the  same  for  two  other  pixel 
thresholds  (70%  and  80%  below  the  maximum  intensity).  Figures  3(d)-3(f)  correspond 
to  hand  measurements  described  earlier.  The  histograms  at  each  axial  location,  obtained 
by  the  two  different  methods,  are  reasonably  close  to  each  other.  To  the  extent  that  the 
histograms  are  not  sensitively  dependent  on  the  axial  location,  we  infer  that  the  buoyancy 
effects  even  at  the  last  station  are  not  severe — even  though  a  priori  considerations  suggest 
otherwise.  This  should  not  come  as  a  total  surprise  because  buoyancy  has  different  degrees 
of  effect  on  different  quantities;  apparently,  the  normalized  histograms  belong  to  the 
category  where  the  effects  of  moderate  levels  of  buoyancy  are  minimal. 

An  attempt  was  made  to  fit  the  histograms  with  a  lognormal  probability  density  function 
using  the  measured  mean  and  standard  deviation.  The  lognormal  fits,  shown  by  the  solid 
curves,  are  reasonable.  The  histograms  were  also  fitted  by  the  Rayleigh  distribution 
function  which  is  similar  in  shape  to  the  lognormal.  The  Rayleigh  distribution  fits  the 
body  of  the  histograms  reasonably  well,  but  not  the  tails.  It  was  found  that  the  lognormal 
estimates  were  closer  to  the  measurements.  Low-order  statistical  data  are  presented  in 
Table  II. 
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Statistical 

averages  front  OH  thickness 

measurements. 

x.  d 

Average 
thickness  (mm) 

rms 

thickness  (mm) 

Skewness 

I  I  at  ness 

38 

4.40 

1.78 

1.26 

5.07 

58 

4.32 

1.79 

1.33 

6.77 

98 

5.22 

3.59 

1.98 

7.1 1 

One  would  like  to  extend  these  statements  to  the  reaction  regions  themselves.  Unfor¬ 
tunately,  as  already  remarked,  because  of  the  slow  recombination  time  and  temperature 
dependence  of  OH  radicals,  the  visualized  OH  regions  are  larger  than  the  reaction 
zone.  The  exact  relationship  between  the  thickness  distribution  of  instantaneous  and 
post  combustion  regions  is  not  evident.  However,  OH  has  previously  been  noted  as  a 
good  diagnostic  of  reaction  zone  structures  (Hanson,  1987;  Vandsburger  et  al„  1988:  and 
Kychacoff  et  al.,  1984).  Experiments  such  as  those  by  Smyth  et  al.  (1985)  on  methane-air 
diffusion  flames  as  well  as  the  study  of  Bilger  (1989b)  and  Keyes  &  Smooke  (1990). 
have  shown  that  peaks  of  OH  concentrations  coincide  with  high  temperature  areas  and 
tend  toward  the  lean  side  of  the  reaction  zone.  Although  our  reservations  noted  above 
remain  in  force,  the  prospect  that  these  results  are  applicable  to  instantaneous  flame 
and  mixing  regions  is  promising. 

5.  CONCLUSIONS 

A  hydrogen  diffusion  flame  was  visualized  by  the  laser  induced  fluorescence  of  OH 
radicals.  Measurements  showed  that  the  thickness  of  OH-containing  regions  is  approxi¬ 
mately  lognormal.  Such  findings  may  be  applicable  to  the  flame  thickness  and  to  mixing 
interfaces. 

ACKNOWLEDGEMENTS 

We  thank  Professor  R.  W.  Bilger  for  his  comments  on  an  earlier  draft.  A  careful  reading  by  a  referee 
improved  the  text.  This  research  has  been  supported  by  the  United  Technologies  Corporation  and  the  Air 
Force  Office  of  Scientific  Research. 

REFERENCES 

Barlow,  R.S.,  Dibble.  R.W.,  Chen.  J.-Y.,  and  Lucht.  R.P.  (1990)  Effect  of  Damkohler  Number  on 
Superequilibrium  of  OF!  Concentration  in  Turbulent  Nonpremixed  Jet  Flames.  Combustion  and  Flume, 
82,  2.15. 

Becker.  FF  A.,  and  Yamazaki,  S.  (1978)  Entrainment.  Momentum  Flux  and  Temperature  in  Vertical  Free 
Turbulent  Diffusion  Flames.  Combustion  and  Flame,  33,  123. 

Bilger.  R.  W.  (1988)  The  Structure  of  Turbulent  Nonpremixed  Flames.  Twenty-Second  Symposium  (International) 
on  CombustionlThe  Combustion  Institute,  p.  475. 

Bilger,  R.  W.  (1989a)  Turbulent  Diffusion  Flames.  Ann  Rev.  Fluid  Mech.  21:  101-135. 

Bilger,  R.  W.  (1989b)  Reaction  Zone  Structure  in  Turbulent  Diffusion  Flames.  Joint  International  Conference 
Australia 'Slew  Zealand  &  Japanese  Sections/The  Combustion  Institute.  University  of  Sydney. 

Hanson.  R.K.  (198b)  Combustion  Diagnostics:  Planar  Imaging  techniques.  Twenty-First  Symposium  (Interna¬ 
tional)  on  CombustioniThe  Combustion  Institute,  p.  1677. 

Keyes,  D.  E..  and  Smooke,  M  D  (1989)  Numerical  Solution  of  2-D  Axisymmctric  Laminar  Diffusion  Flames. 
Combustion  Science  and  Technology.  67,  p.  85. 

Kychakoff.  (}..  Howe,  R.  D..  Hanson.  R.  K.,  Drake,  M.  C.,  Pilz.  R.  W„  Lap.  M  .  and  Penney,  C.  M..  (1984) 
Visualization  of  Turbulent  Flame  Fronts  with  Planar  Laser-Induced  Fluorescence.  Science.  224,  382. 
LaRue.  J  C  .  and  Libby,  R  A.  (1974)  Temperature  and  Intermitteney  in  the  Turbulent  Wake  of  a  Heated 
Cylinder.  Phys.  Fluids,  17,  p.  873. 

Pfefferle.  L  D..  Griffin,  T.  A  .  and  Winter.  M.  (1988)  Planar  laser-induced  fluorescence  of  OH  in  a  chemically 
reacting  boundary  layer.  Applied  Optics.  27,  p.  3197. 


OH  REGIONS  IN  TURBULENT  FLAME 


Seitzman,  J.  M.,  Ungut.  A.,  Paul,  P.  H.,  and  Hanson,  R.  ( 1990)  Imaging  and  Characterization  of  OH  Structures  in 
a  Turbulent  Nonpremixed  Flame,  Twenty-Third  Symposium  I  International )  on  Combustion  The  Combustion 
Institute,  p,  637. 

Smyth,  K.  C.,  Miller,  J.  H  ,  Dorfmann,  R.  C..  Mallard,  W.  G.,  and  Santoro,  R.  J.  (1985)  Soot  Inception  in  a 
Methane/Air  Diffusion  Flame  as  Characterized  by  Detailed  Species  Profiles,  Combustion  and  Flume  62, 
p.  157. 

Sreenivasan,  K.  R„  Ramshankar,  R„  and  Meneveau.  C.  (1989)  Mixing,  Entrainment,  and  the  Fractal 
Dimension  of  Surfaces  in  Turbulent  Flows.  Proc.  R.  Soc.  Lond.  A421,  p  79 

Vandsburger,  U„  Seitzman,  J.  M.  and  Hanson,  R.  K.  (1988)  Visualization  Methods  for  the  Study  of  Unsteady 
Nonpremixed  Jet  Flame  Structure.  Combustion  Science  and  Technolojff,  59,  p.  455. 

Wygnanski,  I.,  Fiedler,  H.  (1969)  Some  Measurements  of  the  Self-Preserving  Jet.  J.  Fluid  Mech  38,  p.  577. 


An  update  on  the  intermittency  exponent  in  turbulence 
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The  issue  of  experimental  determination  of  the  intermittency  exponent  y,  is  revisited  and  it  is 
shown  that  the  “best”  estimate  for  it  is  0.25  ±0.05.  This  "best”  estimate  is  obtained 
from  recent  atmospheric  data,  and  is  based  on  several  different  techniques. 


Ever  since  Kolmogorov1  introduced  the  intermittency 
exponent,  more  than  a  dozen  attempts  have  been  made  to 
obtain  its  numerical  value  from  experiments.  The  resulting 
estimates  have  varied  from  as  low  a  value  as  0. 18  to  as  high 
as  0.7.  (See  Ref.  2  for  a  summary  of  early  measurements, 
Table  IV  of  Ref.  3  for  measurements  until  about  1975  and 
Refs.  4-10  for  experiments  made  thereafter.)  It  is  unsatis¬ 
factory  that  this  exponent,  which  plays  a  moderately  im¬ 
portant  role  in  the  theoretical  framework  of  turbulence, 
should  be  known  with  no  better  certainty.  We  present  in 
this  Brief  Communication  a  brief  critique  of  the  various 
definitions  and  experimental  techniques  used  for  determin¬ 
ing  the  intermittency  exponent,  and  provide  the  best  esti¬ 
mate  for  it  from  recent  atmospheric  data. 

Some  of  this  variability  is  clearly  due  to  differences  in 
the  definitions  of  the  intermittency  exponent,  and  it  is 
therefore  essential  to  describe  all  of  them  briefly.  We  find  it 
convenient  to  adopt  different  symbols  to  denote  the  differ¬ 
ent  definitions  of  the  intermittency  exponent.  First,  there  is 
the  constant  y^  in  Kolmogorov’s  lognormal  hypothesis1 
given  by 

oz=A+yi\og(L/r).  (1) 

Here,  c2  is  the  variance  of  the  logarithm  of  en  which  is  the 
energy  dissipation  averaged  over  an  interval  of  size  r,  L  is 
the  “macro”  or  “integral”  scale  of  turbulence,  and  the  ad¬ 
ditive  constant  A  is  presumed  to  depend  on  the  large  scale 
of  the  flow. 

On  the  basis  of  a  generic  self-similar  cascade  model — 
for  relevant  background,  see  Refs.  2  and  11;  for  a  more 
modern  interpretation  in  terms  of  multifractals,  see  Refs.  9, 
10,  and  12 — we  can  write 

< €?>-/■ (2) 

From  Novikov’s11  argument  (discussed  at  some  length  in 
Ref.  2,  Sec.  25),  one  can  also  write  for  homogeneous  tur¬ 
bulence  that 

<e(*)€(*  +  r)>  (3) 

By  Fourier  transforming  Eq.  (3),  one  obtains  for  the  spec¬ 
tral  density  of  the  energy  dissipation  E\(k)  ~~k]  which 

we  shall  write  as 

(4) 

Finally,  use  has  also  been  made  in  the  past  of  the  definition 
that 

<[fW-<t>][fU+r)-(f)|)~r  "s.  (5) 


Equation  (5)  is  equivalent  to  Eq.  (3)  only  at  very  high 
Reynolds  numbers. 

It  is  safe  to  say  that  — /z3,  and  that  y^=y^  trivially. 
However,  Eq.  (4)  may  not  yield  satisfactory  results  be¬ 
cause  of  crossover  effects  discussed  by  Nelkin;13  these  ef¬ 
fects  are  important  at  all  finite  Reynolds  numbers.  Further, 
the  scaling  region  apparent  in  one  procedure  may  well  be 
different  from  that  in  another.  The  quantity  we  shall  call 
intermittency  exponent  is  y=y2=y-t-  The  constants  and 
y5  are  good  approximations  to  y  only  if  the  Reynolds  num¬ 
ber  is  very  high.  The  constant  y q  is  not  related  to  the  other 
y's  except  through  the  hypothesis  of  some  cascade  model. 

We  are  now  in  a  position  to  comment  on  experiments 
in  more  specific  terms.  Most  of  the  early  estimates  of  the 
intermittency  exponent  were  made  by  using  Eq.  (4). 
Nelkin13  pointed  out  the  difficulty  in  estimating  y  via  p4 
when  the  scaling  range  is  finite  and,  with  this  background 
and  a  model  calculation,  deduced  a  value  of  0.25  from 
Kholmyansky’s14  spectral  data.  Antonia  et  al.s  followed  up 
with  an  immediate  experiment  in  a  fairly  high  Reynolds 
number  jet  (RA  =  966),  and  confirmed  Nelkin’s  estimate. 
They  obtained,  by  means  of  Eq.  (3),  the  exponent  to  be 
around  0.2.  This  latter  estimate  was  confirmed  by  Ansel- 
met  etal,b  who  also  used  Eq.  (3).  It  is  interesting  that 
Anselmet  et  al.  found,  from  the  same  data  at  a  moderate 
Reynolds  number,  that  y5  was  about  0.45,  emphasizing  the 
finite  Reynolds  number  effect  to  which  we  have  already 
aliuded.  the  recent  multifractal  scaling  work4'10  has 
yielded  a  value  of  about  0.25  ±0.05.  These  same  experi¬ 
ments  also  showed  that  the  measured  value  of  y  j  was  about 
0.23,  in  close  agreement  with  0.25.  This  may  indicate  that 
some  kind  of  cascade  process  is  not  a  bad  model  for  tur¬ 
bulence. 

Despite  this  seemingly  conclusive  evidence  that  the  in¬ 
termittency  constant  is  about  0.25  ±0.05,  older  estimates 
have  often  found  favor  in  the  literature.  For  example. 
Baker  and  Gibson,7  and  Gibson8  have  used  a  value  of  0.5 
for  y.  Gibson  (private  communication)  has  indicated  how 
the  smaller  value  of  0.25  would  imply,  when  used  in  con¬ 
junction  with  the  lognormal  model  for  the  energy  dissipa¬ 
tion,  an  astronomical  value  for  the  energy -containing  scale. 
It  seemed  that  an  effort  at  revisiting  the  issue,  and  settling 
it  conclusively  if  possible,  was  worthwhile. 

We  systematically  examined  the  various  methods  of 
determining  the  intermittency  exponent  using  some  recent 
data  obtained  in  the  atmosphere.  Measurements  were 
made  in  the  atmospheric  surface  layer  about  6  m  above  a 
long  stretch  of  a  wheat  field  canopy,  and  at  a  height  of  2  m 
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FIG.  t.  The  log-log  plot  of  <e(.x)e(x+r) )  vs  r  (assuming  Taylor's  hy¬ 
pothesis).  In  this  and  the  following  two  figures,  r  is  expressed  in  terms  of 
sampling  units=  1/6000  sec  and  the  wind  velocity  was  6  m/sec.  The  best 
estimate  of  n=0.25.  This  so-called  best  estimate  was  obtained  by  trying 
various  values  of  the  exponent  ft,  and  deciding  upon  the  value  of  ft  which 
yielded  the  best  flat  line  in  plots  of  (ordinate/r  '')  in  Eq.  (3). 


above  the  roof  of  a  four-story  building.  The  mean  wind 
velocity  varied  between  4  and  6  m/sec.  Velocity  fluctua¬ 
tions  were  measured  using  a  standard  hot  wire  ( 5  ^m  in 
diameter,  0.6  mm  long)  operated  on  the  constant  temper¬ 
ature  mode  on  a  DISA  55M01  anemometer.  The  anemom¬ 
eter  voltage  was  digitized  on  a  12-bit  digitizer  at  a  sampling 
frequency  of  6000  Hz.  This  sampling  frequency  was  found 
to  be  adequate  for  capturing  most  of  the  small-scale  fluc¬ 
tuations.  The  linearization  of  the  signal  yielded  a  time  trace 
of  the  streamwise  velocity  fluctuation  «(f ).  The  dissipation 
field  was  approximated  by  ( du/dt )2.  This  assumes  that  the 
space  derivative  could  be  approximated  by  the  time  deriv¬ 
ative  according  to  Taylor’s  frozen  flow  hypothesis,  and 
that  one  component  of  dissipation  is  an  adequate  represen¬ 
tation  of  the  total  dissipation  statistically.  It  is  believed  that 
these  approximations  were  not  critical  to  the  determination 
of  the  intermittency  exponent.  The  internal  Reynolds  num¬ 
ber  based  on  the  root-mean-square  velocity  and  the  Taylor 
microscale  varied  between  1500  and  2000. 

The  following  characteristics  of  the  data  reassured  us 
about  their  quality.  The  power  spectral  density  showed  a 
sizeable  scaling  range  with  a  slope  of  approximately  —  5/3 
(actually  slightly  steeper).  The  third-order  structure  func¬ 
tion  varied  linearly  with  the  separation  distance15  for  more 
than  a  decade  of  variation  in  the  separation  distance. 

The  most  unambiguous  methods  of  determining  p  are 
from  the  Eqs.  (2)  and  (3);  to  our  knowledge,  all  previous 
measurements  that  used  these  methods  have  given  a  value 
close  to  0.25.  Figures  1  and  2  from  the  atmospheric  data 
show  that  both  Eqs.  (3)  and  (2)  yield  essentially  the  same 
answer,  namely,  p  is  about  0.25.  (Note,  however,  that  the 
scaling  ranges  in  the  two  figures  are  somewhat  different.) 

While  the  experience  of  Anselmet  et  al.k  would  suggest 
that  there  is  little  point  in  using  Eq.  (5),  we  obtain  more  or 
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FIG.  2.  The  log-log  plot  of  (e;>  vs  r  for  the  same  set  of  data  used  in  Fig 
1.  The  best  slope  is  0.27. 

less  the  same  exponent  irrespective  of  whether  we  use  Eqs. 
(4)  or  (5).  See  Fig.  3.  We  tentatively  conclude  that  the 
discrepancy  in  the  two  methods  of  measurement  in  the 
experiments  of  Anselmet  et  al.  is  their  lower  Reynolds 
number. 

We  now  use  an  alternative  method  for  the  determina¬ 
tion  of  the  intermittency  exponent.  This  method  uses  the 
so-called  multiplier  distributions4111-  appropriate  to  the 
inertial  range.  Briefly,  the  multipliers  M a  are  defined  as  the 
ratio  of  the  energy  flux  (or  the  energy  dissipation  rate) 
contained  in  an  eddy  of  size  r  in  the  inertial  range  to  that 
in  a  box  of  size  ar,  where  a  is  an  integer  that  equals  the 
number  of  subeddies  into  which  a  parent  eddy  is  presumed 
to  break  up.  The  multipliers  Mu  are  random  variables  that 
possess  a  well-defined  probability  density  function,  p{Ma). 
In  Ref.  12,  it  was  shown  that,  for  any  given  a,  p(Ma)  is 
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FIG.  3.  <[<?(*)  <r)][f(x  f  r)  <()]>  vs  r  for  the  same  set  of  data  used 

in  Fig.  I.  The  best  slope  is  0.27. 
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independent  of  r  for  all  r  in  the  inertial  range.  Chhabra  and 
Sreenivasan1"  obtained  the  distributions  of  these  multipli¬ 
ers  for  several  values  of  a,  and  published  results  for  a-  2, 
3.  and  5  (Fig.  2  of  Ref.  12).  The  fact  that  p{Ma)  is  inde¬ 
pendent  of  r  for  any  given  base  a  can  be  used  to  advantage 
in  obtaining  converged  values  of  high-order  moments  that 
one  would  otherwise  not  be  able  to  measure  reliably.  Al¬ 
though  the  probability  density  p(Mj  depends  on  a,  the 
various  scaling  exponents  given  by  the  multifractal  spec¬ 
trum  are  independent  of  a.  In  particular,  the  intermittency 
exponent  is  independent  of  the  precise  value  of  a  chosen. 

It  is  easy  to  compute  p2  from  multiplier  distributions. 
Let  us  start  with  an  eddy  of  size  unity,  and  normalize  the 
total  energy  dissipation  in  that  eddy  to  be  unity.  Let  such 
an  eddy  break  down  into  a  pieces,  and  let  the  energy  dis¬ 
sipation  contained  in  each  piece  (as  a  fraction  of  the  total 
contained  in  the  parent  eddy)  be  picked  randomly  from 
the  appropriate  distribution  p(Ma).  Let  us  examine  the 
subeddies  that  have  resulted  after  n  steps  in  this  cascade. 
By  the  definition  of  the  multipliers,  the  amount  of  energy 
dissipation  Er  contained  in  a  box  of  size  r—a~  "  is  given  by 


Er= 


i=n 


n  Mau). 

/= t 


(6) 


Here,  by  the  respective  definitions  of  the  two  quantities, 
Er=r  er  (in  one  dimension).  The  product  in  Eq.  (6)  over 
the  index  /  (which  varies  between  1  and  the  current  step  n 
of  the  cascade)  occurs  because  the  total  energy  dissipation 
E,  contained  in  a  subeddy  of  size  r—a  n  is  the  product  of 
the  multipliers  from  one  step  to  another  starting  with  the 
scale  unity.  These  multipliers  will  be  chosen  from  the  ap¬ 
propriate  probability  density  function  p(MJ. 

It  follows  that 


<(/v)2>=(  n  Ma(n  )  -  n  (md,  (7) 


where  the  products  are  taken  over  all  values  between  1  and 
n ,  and  the  angle  brackets  denote  probability  averages  over 
the  appropriate  multiplier  distribution.  The  last  step  fol¬ 
lows  because  the  multipliers  are  statistically  independent  of 
the  step  of  the  cascade.  Noting  that  Er~r  er=a~"€r  and 
that  II  (Af2)  =  ( (M2a) )",  we  can  show  by  a  few  steps  in¬ 
volving  only  algebra  that 

p=\oga(Ml)  +  2.  (8) 


From  the  probability  distribution  p(Ma)  measured  in  Ref. 
12,  we  obtain  that  /z=0.23±0.05  for  all  three  values  of  a 
(2,  3,  and  5).  This  is  consistent  with  the  present  estimates 
obtained  by  other  methods.  It  should  be  noted  that  similar 
calculations  made  in  Ref.  4  yielded  values  between  0,18 
and  0.22.  Van  Atta  and  Yeh4  also  note  that  these  values  are 
consistent  with  Kholmyanski’s  result  of  0.23. 

Using  the  triangular  approximation  to  p(Ma)  recom¬ 
mended  in  Ref.  12  for  u  =  2,  we  can  analytically  show  that 


(:W~>  7/24,  which  gives  /i  0.22.  For  (hep  model  of  Ref. 

16.  also  for  a-- 2.  (Ml,)  0.29,  which  gives  p  0  2] 

All  these  considerations  lead  us  to  conclude  that  the 
"best  estimate"  for  the  intermittency  exponent  is  ft  0,25 
±0.05. 

Incidently,  all  our  estimates'1  for  the  intermittency 
exponent  in  the  case  of  the  scalar  dissipation  are  m  the 
vicinity  of  0.35. 
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